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Abstract 

Flat directions are a generic feature of the scalar potential in supersymmetric gauge 
field theories. They can arise, for example, from D-terms associated with an extra 
abelian gauge symmetry. Even when supersymmetry is broken softly there often 
remain directions in the scalar field space along which the potential is almost flat. 
Upon breaking a gauge symmetry along one of these almost flat directions, cosmic 
strings may form. Relative to the standard cosmic string picture based on the abelian 
Higgs model, these flat-direction cosmic strings have the extreme Type-I properties 
of a thin gauge core surrounded by a much wider scalar field profile. We perform a 
comprehensive study of the microscopic, macroscopic, and observational characteristics 
of this class of strings. We find many differences from the standard string scenario, 
including stable higher winding mode strings, the dynamical formation of higher mode 
strings from lower ones, and a resultant multi-tension scaling string network in the early 
universe. These strings are only moderately constrained by current observations, and 
their gravitational wave signatures may be detectable at future gravity wave detectors. 
Furthermore, there is the interesting but speculative prospect that the decays of cosmic 
string loops in the early universe could be a source of ultra-high energy cosmic rays 
or non-thermal dark matter. We also compare the observational signatures of flat- 
direction cosmic strings with those of ordinary cosmic strings as well as (p, q) cosmic 
strings motivated by superstring theory. 
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1 Introduction 



Cosmic strings are one- dimensional topological defects that can be formed in the early 
universe [H[2],[3]. They are created if there is a phase transition in which a U(l) subgroup of 
a continuous symmetry is broken. Cosmic strings are stable because they carry a conserved 
topological charge. This charge is integer-valued, corresponding to Hi(U(l)) = Z, and is 
related to the number of times the phase of the U(l) breaking field winds at spatial infinity [I]. 

Unlike other types of topological defects, such as monopoles and domain walls, cosmic 
strings can be formed at a wide range of energy scales after inflation without severely 
contradicting the observed cosmology. The generic problem with topological defects is that, 
on account of their stability, they can easily come to dominate the energy density of the 
universe [3j. For cosmic strings there is an important loophole. Topological stability only 
applies to infinitely long strings. Cosmic string loops do not carry a net topological charge, 
and they can decay into particle or gravitational radiation. Such loops are formed when string 
segments intersect and exchange ends, or reconnect (or sometimes called intercommute) . 
This allows a network of long cosmic strings to regulate its energy by chopping itself up into 
loops which radiate away. Indeed, for a wide range of initial string densities, analytic and 
numerical studies find that the competing processes of string stretching (from the cosmic 
expansion) and loop formation come to balance each other out. The network evolves towards 
a universal scaling solution whose properties are almost fully characterized by the cosmic 
string tension [3 [6], [3, El [9] , independent of the initial conditions. 

The vast majority of work on cosmic strings has focused on the abelian Higgs model, 
in which a U(l) gauge symmetry is spontaneously broken by the condensation of a charged 
scalar field. In this model, the vacuum expectation value (VEV) of the complex scalar field 
determines the mass of the gauge field, my, and the physical scalar Higgs field, m s , through 
the relations 



where g is the gauge coupling, A is the scalar quartic self-coupling, and v is the VEV of the 
scalar. The relative size of my and ms determines how the strings interact. For my < ms, 
parallel strings tend to repel at large distances while anti-parallel strings attract [3]. These 
strings are said to be Type-II, in analogy with superconductors. When my > ms, the 
strings attract for any relative orientation, and they are said to be Type-I. The attractive 
force between parallel Type-I strings allows them to form stable higher-winding modes. 

In most field theories, including the abelian Higgs model, the masses my and ms are 
naturally of the same order. Much of the previous work on cosmic strings has therefore 
dealt with Type-II or weakly Type-I strings. In the present work, we will instead investigate 
the behavior of very strongly Type-I cosmic strings, corresponding to my 3> ms- Our 
motivation to consider the extreme Type-I limit comes from supersymmetry [10J. As we will 
show below, there exist supersymmetric field theories in which my 3> ms arises in a natural 
way when a U(l) gauge symmetry is broken along a flat-direction of the scalar potential. 
Supersymmetry is essential because it ensures that quantum corrections do not lift the flat 
direction. 



my ~ gv 



ms — V\v, 




3 



The key ingredients in our construction, supersymmetry and a new U(l) gauge symmetry, 
are each well-motivated in their own right independently of cosmic strings. Low-energy 
supersymmetry is one of the most elegant ways to explain the large hierarchy between the 
electroweak scale and the Planck scale [10]. It can also provide a candidate for the dark 
matter (DM) in the lightest superpartner particle (LSP), and in its minimal form, leads to 
an excellent unification of gauge couplings. Supersymmetry also plays an important role in 
superstring theories of gravity. Additional local U(l) symmetries arise in many models of 
new physics such as grand unified models and D-brane constructions [11]. In supersymmetric 
models, such symmetries can also help to solve the \i problem |12j . 

A common feature of supersymmetric theories is the existence of directions in the scalar 
potential that are almost flat. To be precise, an almost flat direction is one for which the 
curvature of the potential near the minimum is much smaller than the scale of the (symmetry- 
breaking) VEV at that minimum. Typically, these directions in field space are completely flat 
at tree-level, when only renormalizable operators are included in the potential, but they are 
lifted by higher- dimensional operators, quantum effects, and supersymmetry breaking. As 
long as the supersymmetry breaking effects are both soft and small, the residual approximate 
supersymmetry prevents quantum corrections from destroying the flatness of the potential. 
When a U(l) gauge symmetry is broken along an almost-flat direction, the scalar excitation 
around the VEV along the flat direction is much lighter than the corresponding massive 
gauge boson. We will show that the cosmic strings associated with this pattern of gauge 
symmetry breaking are of the strongly Type-I sort [131 HH US EH El EH1 ISDjEl 

The interactions and cosmological consequences of strongly Type-I strings can be quali- 
tatively different from those of Type-II and weakly Type-I strings [251 [26J, [27] . When Type-I 
or Type-II cosmic strings intersect, they can reconnect or pass through each other. There is a 
third possible outcome when a pair of strongly Type-I strings intersect. Due to their mutual 
attraction, two strong Type-I strings with topological charges N\ and N 2 can combine to 
form a single stable string with topological charge N zip = (JVi + N 2 ) or N zip = \Ni— N 2 \. At 
the point of intersection, the incident strings can coalesce into a single higher-winding string, 
which may then proceed to grow like a zipper [26J. If this growth continues indefinitely, the 
outcome will be a single higher-winding mode string of horizon length. For Type-II and 
weakly Type-I strings, previous calculations and simulations predict that the outcome of a 
string intersection is reconnection with a probability close to unity, P r ~ 1 [281 [2HJ HES [31]. 
Since reconnection is essential to the formation of string loops, which in turn are essential for 
the strings to be cosmologically viable, deviations away from P r ~ 1 can significantly alter the 
picture of cosmic strings in the early universe. In particular, if string zippering is common, 
there can exist a stable population of higher winding mode strings as well [321 [331 EH [35] . 

Many of the exotic properties exhibited by the strongly Type-I cosmic strings arising 
from supersymmetric flat directions are also found in the (p, q) cosmic strings emerging 
from superstring theory [361 EH EH EHl EH HH H21 H3], consisting of p fundamental F- 
strings and q D-strings. These cosmic superstrings can merge to form the equivalent of 

1 Let us also emphasize that the cosmic strings arising in general (approximately) supersymmetric theories 
need not be associated with a flat direction, and can also be of the Type-II variety. For examples, see 

Refs. miiii nam]. 



4 



higher winding modes. In many cases they also have reconnection probabilities much less 
than unity, P r < 1. However, flat-direction strings differ greatly from these (p, q) strings in 
their microscopic properties. This is borne out in the the relationship between the (effective) 
topological charge and the string tension, as well as in the selection rules for string zippering. 
It may therefore be possible to distinguish (p, q) strings from flat-direction strings with the 
observation of several string lensing events, each with a different apparent relative value for 
the string tension. 

In the present work we study the properties and implications of cosmic strings derived 
from the breakdown of a U(l) gauge symmetry along a supersymmetric flat direction. We 
begin in Section [2] by studying the internal structure of flat-direction strings. Here, we 
present a simple toy model for the flat direction breaking, and we investigate approximate 
solutions to the equations of motion and study the string tensions using variational methods. 
In Section [3] we discuss the interactions between cosmic strings. We apply these results in 
Section HJ where we study the formation and evolution of flat-direction string networks in 
the early universe. The observational signatures produced by these networks will be the 
subject of Section [5j Finally, Section [6] is reserved for our conclusions. 

Several earlier papers have investigated cosmic strings associated with flat-directions [T3J 
[III EE [ISJ 13 EE [2D]. These studies have predominantly focused on the lowest (N = 1) 
winding mode. We expand on these studies by exhibiting an explicit and natural field theory 
model for the strings, and by discussing the new features that arise from the existence of 
stable higher (N > 1) winding modes. These modes significantly alter the cosmological 
picture of the strings. 

2 String Profiles and Tensions 

To begin, we introduce a simple class of models for a supersymmetric flat direction that 
could arise if there exists a U(l) gauge group in addition to those contained in the minimal 
supersymmetric standard model (MSSM). Within these models, we study the cosmic string 
solutions they support. In particular, we find approximate solutions to the classical equations 
of motion subject to the boundary conditions appropriate to a cosmic string, and we use 
these solutions to motivate a variational estimate of the string tension. Even though we focus 
on a particular class of models in the present section, we expect that many of the qualitative 
features that we find are also applicable to other cosmic string solutions associated with flat 
directions. 

2.1 (a, b) Flat Directions 

As a prototypical model for U(l) symmetry breaking along a supersymmetric flat direction, 
we consider the (a, b) model discussed in Ref. [44J. The model consists of a supersymmetric 
£7(1) gauge theory containing chiral superfields $ a and with integer charges a and —b 
respectively. Except for the special case a = b = 1 [45], we will assume that a and b are 
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relatively prime with a + b > 3. Aside from the (1, 1) model, other charged fields must 
be present in the theory for anomaly cancellation. However, these will decouple from the 
present discussion as long as they do not develop VEVs. 

When the charges a and b are relatively prime, the leading superpotential operator built 
from <3> a and is 

where M* is a large mass scale above which our effective theory breaks down. We also include 
the soft supersymmetry breaking operators 

3 -m 2 a \ip a \ 2 - m 2 b \ip. b \ 2 - f -A— tp^ + h.cA , (2.2) 

where ip a and ip-b are the scalar component fields of $ a and <&_&, and A is a dimension-one 
coupling on the order of the soft supersymmetry breaking scale, A ~ \/| m al ~ \/l m bl@ 111 
writing this expression, we have implicitly redefined the scalar components of $ a and 
such that A is real and positive. We have also taken the soft masses for (p a and to be 
tachyonic. 

The leading contributions to the scalar potential in the model are therefore 

Vf = J^reO^rV^r + ia^^ 1 ! 2 ), (2-3) 



i. 

2 



Vd = 9 4 {a\Va\ 2 ~ b\ V - b \ 2 ) 2 , (2.4) 



V sof t = -m 2 a \ip a \ 2 - m 2 b \tp- b \ 2 - { -A^ vlv'-b + h - c ^j 



(2.5) 



With A real and positive, there will be a global minimum of the potential with both (p a and 
ip_b real and positive. This minimum is unique up to gauge rotations. 

If a + b > 3 the potential will be almost flat along the D-flat direction defined by 

a\va\ 2 = b\v-b\ 2 - (2.6) 

Along this direction, the potential is destabilized at the origin, and is only restabilized at 
large field values by the higher dimensional F term operators. Near the minimum, the 
excitation along the flat direction is much lighter than the excitations orthogonal to it as 
well as the gauge bosons. This allows us to integrate out the heavy modes and obtain an 
effective potential for the light excitation. 

Let us restrict ourselves to the flat direction by setting 

ip a = v cos a, ip_b = v sin a, (2.7) 



2 A simple spurion analysis indicates that other, non-holomorphic supersymmetry breaking terms from 
insertions in the Kahler potential are subleading [33]. 
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where 



cos a 



b 

a + 6' 



sin a 



a + b 



The scalar potential for v becomes 

2Q 



V(v) 



with 



-Pv z - 
P = 

Q = 

R = 



a + b 



v 



2\(a+b)/2 



R 



a + b- 1 



2\a+b-l 



(2i 



(2.9) 



bm 2 a + am\ 



a + b 
A 



Mi 



a+b-3 



alb 



a a b 



_(a + b) 



a+b-2 



1/2 



(2.10) 



IAI 



M 



2a+2fe-6 



>b b 



(a + b) a 



+b~2 



In terms of these variables, the minimum is given by 



^ (Q + V / Q 2 + 4Pi? 



(a + b- 1). 



1/ (a+b-2) 



Parametrically, this is on the order of 

«~( ffl Mf- 3 ) 1/( " +M 



(2.11) 



(2.12) 



where m is the generic soft mass. Thus, we expect m<Cu< M*. The true minimum of the 
potential does not lie precisely along the flat direction if m\ ^ m^. However, the deviation 
is very small, and can be expanded in powers of (m 2 a — ml)/g 2 v 2 1. We will discuss this 
further below. 

For the special (1,1) case with field charges ±1, we disallow the bilinear term as in 
Ref. jl5] and only include the next-to-leading order term in the superpotential, 



A 



The various terms in the potential are therefore 



Vf 
V D 

V so ft 



4IA 



Af2 (l^-il 2 + l^-i| 2 ) 

9 (\ |2 I |2\2 



-m\\(fx 



2 2 I I 

- m_ x \ip_i\ 



(2.13) 

(2.14) 

(2.15) 
(2.16) 



In the following sections we will analyze in detail the equations of motion resulting from this 
scenario. 
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2.2 Equations of Motion and Approximate Solutions 



The equations of motion for the system are 

dV 

= + — , (2.17) 

= duF^ + igY^QiitfDpipi), (2.18) 

i 

where = + i g Q is the gauge-covariant derivative. 

To obtain an approximate solution to these equations that describes a cosmic string, it is 
convenient to introduce an Ansatz for the vector and scalar fields. Our Ansatz for a string 
with winding number N is 

<p a = v (1 + e) cos a e iNa * f a (r), 
<p- b = v(l-e)smae- mb *f b (r), (2.19) 

N ~, \ 
As = — a(r). 

gr 

In these expressions, r and are the radial and angular cylindrical coordinates relative to the 
string axis, v is the vacuum expectation value, and cos a and sin a are defined in Eq. (12. 5j) . 
The dimensionless parameter e characterizes the deviation from D-flatness at the absolute 
minimum, and will be treated as a small number. The functions f a (f), /f>( r ), and a(r) are 
undetermined string profiles. They are subject to the boundary conditions 

fa, h, a -> as r -> 0, f a , f b , a -»• 1 as r -> oo. (2.20) 

The relative winding numbers of ip a and allow for both D$ip a and D^_ b to fall off more 
quickly than 1/r as r — >• oo. This is a necessary condition for the string tension to be finite. 

Inserting the profile functions into the equations of motion, we obtain 

= /: + V a - ^ (1 - 5) U - a (-^r) [(1 + eff a - (1 - effl] f a (2.21) 
r r z \a + o J 

1 1 c -iNa4,9V 



v(l + e) c a g 2 v 2 d(fl 
f'l + \f'b - ^ (1 - S) h + b (^) [(1 + eff a - (1 - e ) 2 / 2 ] h (2.22) 



v(l — e) s a g 2 v 2 d(p*_ t 
= a" - la' + (^) [a (1 + e) 2 f 2 + 6(1- e) 2 / 2 ] (1-5). (2.23) 

In these expressions we have separated out the D-term part of the potential by defining 
V = (V — Vd)- We have also written the cylindrical radial coordinate r in units of (gv)^ 1 . 
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Thus, when we discuss the value of r in absolute terms, it will always be relative to the 
scale (gv)~ x . The equations of motion are complicated and non-linear, but we can obtain 
approximate solutions in the three regions r < 1, 1 < r < ju/m, and r ^> gv/m. We 
consider each of these regions in turn. 

Region I: r<Cl 

For r C 1, we expect f a , fb, and a to all be small. Expanding the equations of motion 
to linear order in the profiles, we find 



fa 
fb 



,\Na\ 

i 

,\Nb\ 



(r < 1) 



a ~ r 



(2.24) 
(2.25) 
(2.26) 



This behavior agrees with the expectation from Refs. [2} [3] 
Region II: 1 <C r <C gv/m 

In the intermediate region 1 <C r <C gv/m, we expect f a , fb, and a to all be on the 
order of unity. As we will discuss below, in this region it is also self-consistent to neglect the 
contribution of V = (V — Vjj) to the equation of motion and to set e = 0. The equations of 
motion for f a and fb simplify if we rewrite them in terms of f+(r) and f-{r), defined by 



/+ 
f- 



Kfa + fb) 
(fa — fb) 



fa 
fb 



U + §/- 

/+ - \f- 



(2.27) 



The equations of motion for f a and fb then imply 



o ~ /* + -/!- 

r 

o ^ fl + -f' + -\ 

r 2 



2ab 
a + b 
2ab 



a + b 



(a + 6)/ + + -(a-6)/_ 
2-b)f + + ~(a + b)f„ 



As r grows larger than unity, the boundary conditions imply / + 
slowly varying in this region, the approximate solution for /_ is 



/+/-, 

/+/-, 
1 and /_ 



(2.28) 
(2.29) 
0. If /+ is 



/_ ~ K {V2abf + r) 



-V2abf+r 



(2.30) 



Thus, /_ falls off quickly, corresponding to the damping of the scalar excitation orthogonal 
to the flat direction. With /_ very small, the equation for / + reduces to 



The corresponding solution is 







/+ + "/+■ 

r 



(2.31) 



/oln - 



(2.32) 
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for some constants fo and r$. Our approximate result is self-consistent because /+ is indeed 
a slowly- varying function of r. 

We can also use this expression for /+ to check the range of r over which we can safely 
neglect the effects of the V term in the equation of motion. For /_ < 0, / + ~ 1, this term 
is on the order of (m 2 /g 2 v 2 ) /+, where m is the scale of the soft super symmetry breaking 
terms. The necessary condition for ignoring the V contribution to the equation of motion to 
the level of approximation we are working to is 

SI, -f + »(^)h =* r<C»-l. (2.33) 

To track the evolution of the gauge profile it helps to define 5a — 1 — a. The corresponding 
equation of motion is 

~ 5a" - ^5a - ab fiSa, (2.34) 

where we have made use of the fact that /_ is expected to damp out quickly and that f C 1. 
The solution is 

5a oc r K 1 (V2a~bf+r^ ~ ^ e~^ h ^ r . (2.35) 

Therefore, 5a is damped out exponentially as well, and a quickly approaches unity. Let us 
point out that the physical gauge boson mass is \Jlabgv. Thus, this mass controls the width 
of the gauge field profile (remembering that r is expressed in units of 1/gv here), as well as 
the width of the profile of f-(r). 

Region III: r 3> gv/m 

In the very large field region, r ^> gv/m, the flat potential V and the deviation of e 
from zero become relevant to the evolution of / + and /_. For these large values of r, it is 
convenient to write 

5f + = l- f + , (2.36) 

since we expect \5f+\ 1. Consider first the effect of V and e ^ on the evolution of /_. 
The equation of motion to linear order in /_ and 5f + becomes 

= fl + -f. - [2ab + 0{e)\ /_ - 4 ab e + ( ^ ~ f b ). (2.37) 
r \ g z v z J 

To be able to impose /_ — > 0, we must choose 

1 / 'm 2 . — ml\ 
e =4^(^V^j- < 238 » 

This is consistent with our previous assumption that e C 1. 
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Inserting this value of e into the linearized equation of motion for 5f + , we find 

= dfl + Uf+ - m 2 s Sf + , (2.39) 

where m| is a positive constant on the order of m 2 /g 2 v 2 . In the units we are using, this 
is of the same size as the mass of the light excitation about the almost-flat direction. A 
possible constant term in Eq. ( 12.391) vanishes through the minimization condition for v given 
in Eq. ( 12. lip . The solution for 8f+ in the very large r region is therefore 



5f+ cc K (m s r) ~ / e ~ m sr . (2.40) 
w 2m s r 



Again, this is consistent with the results of Refs. [21 E]- 



2.3 String Tensions 



Having obtained approximate expressions for the string profiles, we estimate the tension of 
cosmic strings in the (a, b) model for various values of the winding number N. Using the 
Ansatz of Eq. (12.191) . the contributions to the tension of a string in the (a, b) model are 



AW/ 7T V 2 



drr 



f^mag / 



IT V 



2N 2 ab / dr 



a + b 
1 



U'aY + 



a + b 



i\2 



(fb) 



(2.41) 



a + b 



a + b 



fi 



7T V 



N 



z>\2 



Upot/n V 2 



dr — (a 



drr -^-J v (fa,fb)- 



Except near the origin, and certainly whenever the potential is relevant, it is a very good 
approximation to set f a — fb — /+■ In this limit, the potential can be written in the form 



g 2 v 2 



V(f) * Si (fl - 1) 



S2 V 



a+b 



D + 



5 2 +6i 



(./- 



2a+2b-2 



1). (2.42) 



Here, we have implicitly assumed that a + b > 4. The dimensionless constants 5\ and 62 are 
given by 



Si 



a + b 
1 



bra 2 + am? 



g 2 v 2 



(2.43) 



A 



aub 



a a b 



[{a + b) 



a+b-2 



1/2 



2+fe-2 



g2 v 2 M a+b-3 

Using the parametric value of the VEV given in Eq. (I2.12p . these constants are of size 

2(o+5-3)/ (a+b-2) 



Si: 



m 
Ml 



(2.44) 
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For M* ~ Mpi and m ~ TeV, we find 10 -30 < 8\ i2 < 10~ 15 . Although the expressions 
presented above were formulated for strings in the (a, b) theory, they can also be applied to 
(1, 1) theory cosmic strings. The correct formulae for the (1, 1) case are obtained by setting 
a = b = 1 in the radial and angular components of the tension (fi ra d and fi ang in Eq. (I2.4ip ). 
but a + b = 4 in the expression for the potential (Eq. (12.421) ). This adjustment accounts for 
our inclusion of terms beyond the leading order for (1, 1) strings. 

To estimate the string tensions, we have used variational methods as in Ref. [13]. Our 
trial profile functions are inspired by the approximate solutions found above. They are 

r < ri 
ri < r < r 2 
r > r 2 
r < r\ 

r 1 < r < r 2 ( 2 -45) 
r > r 2 

r <r a 
r > r a 

The undetermined parameters are {n, r 2 , r 3 , r a , pi, p 2 , P3, p^, p{\. We fix four of them, pi, 
P2, P3, and £>4, by requiring continuity at r = r\ and r 2 , and differentiability at r 2 (where the 
solution is expected to be slowly varying) but not at r%. 

For a (1, 1) model string with winding number N = 1 and 8% = 5 2 = 1 x 10~ 20 , our 
variational estimate of the tension is 

Vrad/nv 2 = 0.09093, 

/W 7 ™ 2 = 0.00247, 

Vmag/KV 2 = 0.00228, (2.46) 

fi pot /nv 2 = 0.00228, 

fi tot /7cv 2 = 0.09796. 

The corresponding values of the variational parameters are 

ri = 14.01, 

r 2 ~ r 3 = 3.112 x 10 9 , (2.47) 
r a = 36.26, 
p 5 = 0.04713. 

Recall that we express all dimensionful quantities in units of 1/gv. 

As expected, the gauge profile is much narrower than the scalar profiles (i.e., r a <ti r 2 ), 
which have substantial support out to r ~ 1/ y/5\^- The small r power-law form of the scalar 
profiles extends out about as far as the gauge profile (i.e., r\ ~ r a ), after which it continues 



{p 1 (r/r 1 y Na \ 
P5+P3 In (^-J 

r p 2 (r/r 1 )\ m \ 
fb(r) = I p 5 +P3 In ^) 
[ 1 -p 4 e- (r ~ r2)/r3 
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Figure 1: Tensions of N — 1 strings as a function of the potential parameters 8\ = 82 for 
various (a, b) theories. 

to grow logarithmically slowly until the profile reaches unity. We also find that the total 
string tension is dominated by the radial contribution. 

To a very good approximation, the shape of the profiles and the value of the string tension 
do not depend on Si and 82 independently, but rather on the combination 

A = 81 + 82/2. (2.48) 

This can be seen explicitly by evaluating n pot using the Ansatz profiles of Eq. (12.451) and 
keeping only the leading terms in the expansion in 1/ ln(5x,2) 1. 

We have investigated a number of other sets of profile functions as well. As long as the 
trial scalar profile increases sufficiently (logarithmically) slowly in the region 1 <C r <C 1/ V8 
and drops rapidly for larger r, we find that the resulting estimates for the string tension are 
very similar. This gives us confidence that our estimates are close to the exact values. 

Fig. [1] shows the dependence of the string tension in the (1,1), (1,3), (2,3), (2,5), and 
(3,4) models on the value of 8\ = 82 = 2 A/3 for a winding number N = 1. Even for very 
small values of £^2, corresponding to extremely flat potentials, the string tension is within 
about an order of magnitude of v 2 . Thus, while the string is very wide in units of 1/gv, the 
VEV still sets the size of the tension. The tensions are also very similar for different values of 
(a, b). This is not very surprising given that the radial portion of the string tension appears 
to be the dominant one. In the i , >1 region, we expect f a ~ fb so that the expression for 
the radial contribution to the tension in Eq. (I2.4ip does not depend explicitly on (a, b). The 
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Figure 2: String tensions as a function of the winding number N for the potential parameters 
Si = 62 = 1 x 10 -20 in various (a, b) theories. Note that the tension of the N = 2 string 
is much smaller than twice the tension of the N = 1 string, thereby allowing stable N = 2 
strings. 

dependence on (a, b) only then comes about through the size of the terms in the potential. In 
generating Fig. [fl we neglected this dependence by specifying the value of Si = S 2 explicitly. 
Our results also suggest that the detailed form of the (non-D) potential does not play a 
significant role in determining the string tension or the string profiles other than to set the 
scale at which the scalar profiles are cut off. 

In Fig. [2] we illustrate the variation of the tension for strings in the (1, 1), (1,3), (2,3), 
(2,5), and (3,4) models with the winding number iV for Si = 82 — 1 x 1CT 20 . These 
tensions increase very slowly with N, approximately logarithmically. As the winding number 
increases, the widths of the vector field profile and the inner portion of the scalar profile do 
too. This allows the angular and magnetic contributions to the string tension to increase 
much more slowly than N 2 . The increase of the profile radii j~i and r a with the winding 
number N is shown in Fig. [3] for a (1, 1) model string with Si = S2 = 1 x 10~ 20 . For both 
ri and r a , the increase with N is very close to linear. The corresponding plots for the other 
values of (a, b) discussed above are nearly identical. Unlike r\ and r a , varying N has very 
little effect on T2- 

We can combine the results presented above into a simple approximate parametrization 
of the string tensions. The string tension increases close to logarithmically with the winding 
number N, but has a more complicated dependence on Si and S2 , primarily through the 
combination A = Si + S 2 /2. In the range 1 < iV < 100, 10 -30 < A < 10~ 15 , and Si and S 2 
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Figure 3: Dependence of the inner scalar profile width (n) and vector profile width (r a ) on 
the winding number number N for a (1, 1) model string with 5i = 62 = 1 X 1CT 20 . 

within an order of magnitude from each other, the tension of a (1, 1) string is reproduced to 
an accuracy of a few percent by the empirical formula 



4.2 14 
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ln(l/A) ln 2 (l/A) 



1 + 



2.6 57 



ln(l/A) ln 2 (l/A) 



InJV . (2.49) 



Since the tension of an (a, b) theory string is very similar to that of a (1, 1) theory string for 
a given set of values of 5± and 62, this formula also provides a reasonable approximation to 
the tension of strings in these more general theories. 

In summary, we find that the cosmic strings that arise from breaking a U(l) gauge 
symmetry along an almost flat direction within the (a, b) models are very strongly of the 
Type-I variety. The qualitative features of these strings can be characterized by two scales: 
the VEV v; and the scale of the curvature near the minimum m, which in the present case 
is set by the soft supersymmetry breaking scale m ~ ^/\m a \ 2 ~ a/I 777 *! 2 ~ A. It is the 
hierarchy m<Cw that makes the potential flat. The tension of flat-direction strings is about 
fi ~ 0.l7ri; 2 , while their total thickness is w ~ m" 1 . The internal structure of the strings 
consists of a thin vector field core, of width close to t> -1 , surrounded by a much broader 
scalar profile of radius m -1 . Flat-direction cosmic strings also have stable higher modes. 
The tension of these modes grows very slowly with the winding number N, increasing as 
IniV with a small coefficient. 

These features are much different from those of ordinary cosmic strings derived from the 
abelian Higgs model, for which the relevant scales are all on the order of the VEV v. On the 
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other hand, the qualitative structure and the tensions of strings derived from the (a, b) model 
presented above are in agreement with other studies of flat- direct ion cosmic strings [T3J [HI 
IT5l [T7] . Within the (a, b) models, we find that the form of the string profile away from 
the central core and the tension can be described well from a knowledge of m and v alone, 
without reference to the precise form of the potential (or a and b). This suggests that many 
of the results of the following sections, where we investigate the phenomenological features 
of (a, 6)-theory flat-direction cosmic strings, will apply to flat-direction strings derived from 
other theories as well. 



3 String Interactions 

When a pair of Type-II or weakly Type-I abelian strings with the same winding number 
intersect, there are effectively two possible outcomes. They can simply pass through each 
other, or they can exchange partners and reconnect (intercommute). When a pair of strongly 
Type-I N = 1 strings collide, there is a third possibility [26] . Studies of Type-I strings in the 
abelian Higgs model suggest that the force between string segments is attractive. Thus, the 
segments can pull together near the intersection point to form a length of N = 2 string, which 
is stable and lower in energy than a pair of N = 1 segments. Under favorable conditions this 
segment will grow, effectively zippering the pair of N = 1 strings into a single N = 2 string. 
When even higher-winding modes of strongly Type-I strings are stable as well, we can also 
consider the outcome of the intersection of two strings with general winding numbers Ni 
and N2. Besides passing through each other, the topology of the configuration permits the 
formation of zippers with winding numbers | N\ + N2 | and | JVj — N2 | ■ 

Reconnection plays an essential role in the cosmological evolution of a cosmic string 
network. It allows the network to modulate its energy by forming string loops, which can 
decay away. Without reconnection and loop formation, the energy density in the string 
network could come to dominant the universe [21 [3]. Analytic estimates and numerical 
simulations of Type-II and weakly Type-I strings in the abelian Higgs model suggest that 
the probability that a pair of strings will reconnect after they intersect is close to one, 
P r ~ 1 [291 EH El]- However, this result need not apply to very strongly Type-I strings. 
These strings can form zippers, and therefore the probability of reconnection in a string 
collision may differ from unity. This can have important consequences for the evolution of a 
string network in the early universe. 

In this section we investigate how flat-direction cosmic strings interact with each other. 
We begin by discussing the forces between a pair of string segments. Next, we study 
the reconnection and zippering of strings when they intersect. Zippering can reduce the 
probability of reconnection, and it can also lead to qualitatively new string structures that 
cannot be formed by Type-II strings. We investigate how these features alter the formation 
of string loops. The results of this section are applied in the sections to follow. 
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3.1 Inter- String Forces 



We found in Section [2] above that the tension of an N = 2 flat-direction cosmic string is 
considerably lower than twice the tension of an N = 1 string. Therefore bringing a pair of 
N = 1 strings together (adiabatically) from infinity to form an N = 2 segment will lower the 
total energy of the system. As a result, we expect the (non-gravitational) force between a 
pair of parallel flat-direction cosmic strings to be attractive. More generally, we expect the 
interstring force to be attractive for any other relative orientation as well. 

Our expectation is supported by both analytic estimates of the interstring forces in the 
abelian Higgs model [1HJ 0H 08 J Jfl as well as in numerical investigations [I9j [50j EE] . It is 
argued in these works that the contributions to the interstring force from the vector profile 
are attractive only for anti-parallel strings and repulsive otherwise, while the scalar profile 
contributions are always attractive. For Type-I strings, the scalar profile is wider than the 
vector profile and its contribution to the force has a longer range and is always dominant. 
The vector profile has a larger range for Type-II strings explaining why the force between 
parallel strings is repulsive. The scalar profile in flat-direction strings is much wider than 
the vector profile, so the results obtained in the abelian Higgs model suggest that the force 
between these strongly Type-I strings is attractive as well. 

An alternative possibility, consistent with the energetics, is that the interstring force 
between flat-direction strings is repulsive at distances larger than the string width, and only 
becomes attractive when the strings overlap significantly. Even if this were true, it would 
likely not have a large effect on how these strings interact in the early universe. Since the 
strings we are studying are local (gauged), the interstring force has a very short range, falling 
off exponentially outside the string core. When a pair of strings approaches an intersection, 
the interstring forces will be non-trivial only in the small region near the intersection point, 
and hence the interaction energy will be finite. We expect the energy required to overcome 
this barrier, if it is present, to be much smaller than the initial kinetic energy carried by the 
incident string segments. 

3.2 String Reconnection and Zippering 

A pair of strings with the same winding number is said to reconnect (or intercommute) if 
they exchange ends upon intersection. The result of this process is illustrated in Fig. 01 
following Ref . [26] . The initial state consists of two infinite straight strings, each with speed 
v and a relative angle a, approaching each other along the z-axis. After exchanging ends, 
causality implies that the segments of the strings very far (spacelike-separated) from the 
intersection point continue along their original trajectories. Connecting these asymptotic 
segments are new segments moving in the ±y directions. The labels 1 and 2 in the figure 
indicate which incident string the corresponding asymptotic string segment came from. The 
total length of string in the final configuration is clearly less than in the initial. Energy is 

3 However, when attempting to reproduce the argument of [47) we found an opposite sign in the scalar 
term at large string separations. Thus, we are not sure that argument is definitive. 
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Figure 4: Pictorial representation of string reconnection in the xy and xz planes following 
Ref. [26J. The initial state consists of string 1 and string 2 approaching each other along 
the z-axis, each with speed v. In the xz plane, we show only the lower string portion. The 
labels 1 and 2 indicate which of the incident strings the corresponding segment was derived 
from. 

conserved because the newly-formed segments carry a velocity v' in the indirections. 

Over distances that are large compared to the string width but small compared to the 
horizon size, the motion of cosmic strings should be well-described by treating them as 
ideal Nambu-Goto (NG) strings propagating in a flat spacetime background. Therefore a 
necessary condition for string reconnection is that the initial and final configurations be 
kinematically allowed in the NG approximation. It is not hard to check that for any initial 
relative velocity v and for any relative angle a (as defined in Fig. HJ), this is the case [26] . 

The existence of a classical string solution for reconnection does not imply that it actually 
occurs whenever a pair of strings intersect. The precise outcome depends on the internal 
structure of the strings, which is highly non-linear and very difficult to treat analytically. 
Much of the work on this topic has therefore consisted of lattice simulations of the cor- 
responding classical field configurations in the abelian Higgs model for Type-II or weakly 
Type-I strings. These simulations generally find that the probability of reconnection in a 
string intersection is close to unity except for very large initial velocities, v > 0.9 [29l [30l [3Tj . 
Early attempts to study this question analytically, by comparing the interaction time of the 
fields in the string core to the time it takes for the pair of strings to pass through each other, 
find much the same result |28j. 

In addition to reconnecting or simply passing through each other, when a pair of strongly 
Type-I strings intersect they can also zipper into a segment with a higher (or lower) winding 
number [25| [26l [27] . This is illustrated in Fig. [5j following Ref. |26j , where the initial state 
consists of two strings with the same winding number ~N\ = N2 = N approaching each other 
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Figure 5: Pictorial representation of string zippering in the xy and xz planes following 
Ref. |26j. The initial state consists of string 1 and string 2 approaching each other along 
the z-axis, each with speed v. In the xz plane, we show only the lower string portion. The 
labels 1 and 2 indicate which of the incident strings the corresponding segment was derived 
from. 

along the z-axis, each with initial speed v. When the strings intersect, a new segment of 
winding number N zip = 2N is formed along the x-axis. This is the zipper. Under favorable 
conditions it proceeds to grow along the x-axis at the speed v zip . The string segments far from 
the intersection point (labelled by 1 and 2 in Fig. [5]) continue along their initial trajectories 
on account of causality. 

Zippering has received much less attention than reconnection, and we know of only a 
handful of simulations that have studied it [25| [271 [52] . If string zippering is efficient, it 
will reduce the probability of reconnection. Given the importance of reconnection for the 
evolution of cosmic strings in the early universe, this is a crucial issue to be resolved. 

As for reconnection, a necessary condition for string zippering is that it be classically 
allowed in the NG approximation. Again, this condition is only a necessary one, and the 
existence of a classical zippering solution does not imply that it actually takes place. Classical 
zippering solutions have been constructed in Refs. [26| [53]. Unlike for reconnection, there 
exist significant kinematic constraints on zippering due to energy conservation. For a pair 
of strings with identical winding numbers N, initial speeds z/, and a relative angle a, the 
kinematic constraint on forming a zipper with N zip = 2 iV is found to be [26] 



where /i^ is the tension of the incident segments and is the tension of the zipper. The 
total length of the zippered configuration is greater than the initial state. Thus, a zipper 
can form only if it tends to lower the energy of the configuration due to the string tension, 
which requires < 2/ijv@ On the other hand, zippering does not occur if the incident 

4 The total energy of the configuration is conserved because parts of the interacting string segments gain 
kinetic energy. 
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strings collide with too great a velocity u, or if the relative opening angle between strings 
with the same winding orientation is too large. 

In Fig. [6] we show the kinematic constraints on the zippering of a pair of N = 1 strings, 
in terms of the incident relative velocity v and the relative angle a, defined in Fig. [5j The 
region where zippering is kinematically allowed lies below the curves. The dashed line for 
weakly Type-I strings was obtained assuming ^2/^1 = 1-9. The solid line corresponding 
to the kinematic constraint on a strongly Type-I flat-direction string was obtained using 
the tensions from Eq. (I2.49p . and found to be ^2/^1 — 1-06. As we will discuss below, the 
typical relative velocity of a pair of strings in the early universe is expected to be less than 
about v < 0.7. Thus, zippering of flat-direction strings in the early universe is kinematically 
allowed for a wide range of relative angles. In the weakly Type-I case, zippering is only 
possible for small relative velocities and angles making it much less likely to occur. This is 
why flat-direction strings can have a qualitatively different behavior in the early universe 
from the strings in the abelian Higgs model. Recall that there are no kinematic constraints 
on reconnection. 




Figure 6: Kinematic constraints on zippering of two N = 1 strings to form an iV = 2 
string, in terms of the angle a indicated in Fig. [5] and the relative velocity v. The allowed 
regions lie below the curves. The dashed line corresponds to weakly Type-I strings, with 
A*2/Vi — 1-9- The solid line corresponds to strongly Type-I strings associated with a flat 
direction potential, with A = 10~ 20 , and tensions computed according to Eq. (12.491) . which 
gives ^2/^1 — 1.06. 

More generally, zippering can occur between Type-I strings with different tensions. Inci- 
dent strings with winding numbers Ni and N 2 can zipper into segments with N zip = (N1+N2) 



20 




Figure 7: Kinematic constraints on zippering of strongly Type-I strings, for some examples 
involving higher winding numbers. The allowed regions lie below the curves. We have taken 
A = 10~ 20 , and tensions computed according to Eq. (12. 49ft . 

or N zip = \Ni — N 2 \ [25J. The corresponding kinematic constraint for the zippering of strings 
with unequal tensions was deduced in Ref . [53] . Zippering is only possible when the tension 
of the zippered segment is less than the sum of the tensions of the incident segments. Even 
when this condition is met, zippering is only allowed for a limited range of relative incident 
velocities v and relative angles a (as defined in Fig. [5]). We illustrate these kinematic 
constraints in Fig. [7] for the incident string pairs N\ = 1 and N 2 = 2, N\ = 1 and N 2 = 100, 
and Ni = 100 and N 2 = 101. The tensions of these strings were computed using Eq. (12.491) 
with A = 10~ 20 , which applies to flat-direction strings in the (a, b) theory described in 
Section [21 As before, the regions in which zippering is kinematically allowed lie below the 
curves. The kinematic constraints on flat-direction strings are not overly restrictive, and 
zippering of various sorts is possible over a wide range of relative velocities v and relative 
angles a (as defined in Fig. [5]). 

When a pair of strings with winding numbers N\ and N 2 intersect, they can pass through 
each other, or they can form a zipper with N Z i P = (JVi + N 2 ) or \N\ — N 2 \. If Ni = N 2 , 
these strings can also reconnect H There is no kinematic restriction on reconnection, and the 
kinematic constraints on zippering (into one of | N\ ± N 2 1 ) are fairly mild. Having determined 
the possible outcomes, it is a much more difficult task to determine which of them actually 

5 In fact, string reconnection can be treated as the formation of a zipper with N zip = 0. The classical NG 
zippering solution reduces to the reconnection solution in this limit. The absence of a kinematic constraint 
on reconnection can be seen by setting \i Z i V = in Eq. (|3 . 1 1) . 
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occurs. The answer depends on complicated non-linear field dynamics within the string 
cores, and would appear to be tractable only through numerical simulation. Unfortunately, 
even this approach is further complicated by the large disparity in scales between the sizes 
of the vector and scalar profiles within the strings. Such a simulation is beyond the scope of 
this paper. 

To proceed, we will assume that zippering or reconnection are likely to occur when 
they are kinematically allowed. Given the high probability of reconnection of abelian Higgs 
strings, this assumption does not seem overly optimistic. When both zippering and recon- 
nection are possible, or when more than one kind of zippering is allowed, we will make use 
of the fact that the net force between a pair of strings is expected to be attractive. This 
suggests that, near the intersection point, the strings will pull together in whichever way is 
easiest. Thus, for a pair of strings with winding numbers Ni and N 2 , we will assume that 
a zipper with N zip = Aq + N 2 forms when a < 45° (provided it is kinematically allowed), 
and that N zip = \Ni — N 2 \ results for a > 45°. We identify the case Aq — N 2 = with 
reconnection. 

Our assumptions are compatible with the two simulations we know of that treat the 
zippering of Type-I (abelian Higgs) strings [251 [27]. In both of these analyses, zippering 
appears to be a generic outcome of a low-speed string intersection. In Ref. [27], the strings are 
found to grow until they reach the size of the box used for the simulation, after which they pull 
apart. This appears to be the result of the boundary conditions applied to the box. We expect 
that in the applications of our assumptions about string zippering and reconnection, our 
qualitative results will still hold true provided the zippering and reconnection probabilities 
are of order unity. 

We end this section with a brief comment of comparison regarding (p, q) cosmic strings 
arising from superstring theory. Like the flat-direction gauge-theory cosmic strings under 
consideration, (p,q) strings are also able to reconnect and form zippers |36j. Even so, there 
are several important differences between the inter-string interactions within these two classes 
of cosmic strings. The reconnection of (p, q) strings is a quantum mechanical process that 
can be related to amplitudes in superstring theory [371 ES EH HOj HH H2J H3] . In this sense, it 
is more tractable than the non-linear classical calculation required for field theory strings. It 
is found that the reconnection probability for (p, q) strings can be much smaller than unity, 
P r ~ 1CT 3 — 1, depending on the underlying microscopic details. The rules for zippering are 
also different for (p, q) strings. An initial state consisting of the modes (p, q) and (p', q') can 
form a zippered state with (|p±p'|, q±q'), which is similar to the topological rule for Type-I 
field theoretic strings presented above. However, a (p, q) cosmic string is stable only if p and 
q are relatively prime integers, and thus the resulting zipper may sometimes decay into lower 
string modes. A recent numerical simulation of a toy-model for (p, q) cosmic superstrings 
has found that long-lived zippered states are formed provided the forces between the strings 
are short-ranged [52] . 
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Figure 8: Two possible ways to form a loop from the self-intersection of a string segment. 
Possibility 1, in which a free loop is formed by string reconnection, can occur for both Type-I 
and Type-II strings. Possibility 2, in which the loop remains connected to the parent string 
by a zippered segment of a higher winding mode string, is only possible for Type-I strings. 

3.3 Loop Formation 

Reconnection plays a crucial role in the evolution of a cosmic string network because it is 
the means by which string loops form. String loops are not topologically stable, and their 
decays transfer energy out of the string network. When cosmic strings are also able to form 
zippers there are new ways for string loops to form and interact. In the present section we 
enumerate some of these additional possibilities. We will discuss the resulting effects on the 
cosmological evolution of a string network in Sections H] and [5j 

In Fig. [8] we illustrate the two ways in which a loop can form when a string intersects 
itself. The first possibility produces a free loop through the reconnection of the intersecting 
segments. This can occur for both Type-I and Type-II strings, and is the standard mechanism 
for loop formation. The loop produced is free from the parent string. The second possibility 
for loop formation through self-intersection involves zippering of the connecting segments. 
The loop formed in this way remains bound to the parent string by a zippered segment of 
winding number N zip — 2Ni, where N\ is the winding number of the parent. We expect 
the zippered segment formed in this way to grow until the opening angle at the junction 
approaches the kinematic bound given in Eq. (13. ip . Subsequently, provided there are no 
disturbances on the string large enough to rip the zipper apart, the bound string loop will 
remain attached to the parent string as it radiates and shrinks to naught. 

String loops can also be formed by the double intersection of a pair of curved strings. 
Suppose the incident strings have winding numbers Ni and N 2 - The topologically-allowed 
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Figure 9: Three ways to form a loop from the overlapping intersection of a pair of Type-I 
cosmic strings with winding numbers Ni and N 2 . In the figure, we have labelled the net 
winding number of each string segment. 

loops that can form in this way are illustrated in Fig. [91 At each intersection, there are 
two ways for the strings to interact with each other by zippering; they can form segments of 
winding number N zip = Ni+N 2 or N zip = Nx — N 2 . (Here and only here, the sign of Ni should 
be understood as specifying the relative orientation of the string segment.) Possibility 1, in 
which both intersections produce segments of winding N zip = iV 1+2 = Ni + N 2 corresponds 
to the usual Type-II outcome when Ni = —N 2 . Possibility 2 has both zippered segments 
with windings N zip = Ni_ 2 = N% — N 2 . It reduces to the standard Type-II case for Ni = N 2 . 
Possibility 3 has zippered segments with winding N± + N 2 and Nx — N 2 . It is not immediately 
obvious how these configurations will evolve, but we speculate that the loops will shrink, 
either through zipper growth or loop radiation, until only a single zippered segment remains. 
The multiple outcomes shown in Fig. also illustrate some of the many new qualitative 
features of a string network consisting of strongly Type-I strings. 



4 Cosmic String Formation and Evolution 

Cosmic strings are much less strongly constrained by cosmology than most other types of 
topological defects (2J [3] . The reason for this is that a network of cosmic strings is able to 
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regulate its energy density by forming loops, which radiate away. Without loop formation, 
the energy density in a cosmic string network would scale as a -2 , redshifting more slowly 
than both matter (a -3 ) or radiation (a -4 ), and could come to dominate the universe. Instead, 
when strings are able to form unstable loops, numerical and analytic simulations suggest that 
the energy density of a string network tracks the dominant background matter or radiation 
density [51 [6j [7J [H [9]. This behavior is called scaling. In the scaling regime, the energy 
density of the string network makes up a fixed proportion of about G\i of the total energy 
density, and this proportion is nearly independent of the initial string density. As long as 
Gfi is not too large, Gfi < 3 x 1CT 7 [63], EU [651 ESI E7J [68], cosmic strings are generally 
consistent with existing cosmological bounds. 

The behavior described above was deduced from the study of Type-II abelian Higgs string 
networks containing only a single string species [21 [3]. Strongly Type-I strings associated 
with super symmetric flat directions can modify this picture in a couple of important ways. 
First, flat-direction strings have stable higher winding modes. Even if modes with N > 1 
are not formed initially, they can be produced as the network evolves by the zippering of 
lower string modes. This opens the possibility that flat-direction strings form a multi-tension 
string network consisting of many different species. The second reason why the evolution of 
flat-direction strings in the early universe is likely to be different than for ordinary strings is 
the flatness of the scalar potential. If the U(l) gauge symmetry corresponding to the strings 
is restored after (or near the end of) primordial inflation, it is likely that there will be a 
second, later period of thermal inflation [5H [56]. Flat-direction strings would be formed at 
the end of thermal inflation, and hence their initial evolution is expected to be significantly 
different from that of abelian Higgs strings. 



4.1 Thermal Inflation and String Formation 

Thermal inflation occurs due to the sensitivity of flat potentials to thermal corrections [HH 
1551 156] . This flatness can be quantified by the large disparity between the size of the 
curvature scale m ~ 10 2-3 GeV and the size of the VEV, v > 10 11 GeV. At the symmetry- 
preserving origin of the field space, there are additional light degrees of freedom. These 
induce significant corrections to the effective potential near the origin, making it stably 
concave at high temperatures, with a curvature scale on the order of the temperature T. For 
m < T < d, a second lower minimum can develop far from the origin, close to the T = 
vacuum. If the system begins in the symmetry preserving phase, thermal corrections will 
trap it at the origin until the temperature falls down to T ~ m [571 ITS"] . While the system 
is trapped at the origin, it has an excess vacuum energy on the order of m 2 v 2 . Once the 
temperature of the universe falls below ^/m v, the false vacuum energy can become dominant 
and drive a period of inflation. 

Thermal inflation lasts only until T falls down to m. The number of e-foldings of 
expansion is therefore less than [56] 

Ni ~ - In (v/m) ~ 10 H — In 



/ v \ /10 3 GeV\ 
V10 14 GeV J V ™> / 
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This is not enough expansion to replace primordial inflation. At the end of thermal inflation 
the system evolves to the true minimum of the potential. In this regime the constituent 
fields $ a and both condense, and the theory can be described in terms of a light 
chiral supermultiplet corresponding to the flat direction along with a heavy massive vector 
supermultiplet |44|. The scalar component of the light chiral multiplet rolls down the 
potential to the true minimum and begins to oscillate. The false vacuum energy is transferred 
to the energy of the oscillations, which redshifts like matter, and dominates until the scalar 
field decays into radiation and reheats the universe. 

The reheating process can be described by the system of Boltzmann equations 

= -ZHp^ - T^p^, (4.2) 
p r = -AH p r +T < j > p (j) , (4.3) 

where p^ is the energy density of the scalar field oscillations, p r is the energy density in 
radiation, is the scalar field decay rate, and the Hubble constant H is given by 



a 811G .. .. 

H=- = \— Ptot . (4.4) 
a V o 

Here, p tot is the total energy density in the universe. During reheating, p tot is dominated by 
p r and p^. The initial values for these evolution equations are p r ~ m 4 , p^ ~ m 2 v 2 , and 
t{ ~ 10 if" 1 ~ lOMpy/mv. The generic value of the flat-direction decay rate is [56] 

TfL 

T^ = 7— , (4.5) 

v 

with 7 a constant less than or of order unity. Once the scalars decay at about the time 
tRH — r^ 1 , the universe becomes radiation dominated with a reheating temperature of 

T RH ~ g-V^Mnr) 1 ' 3 (4.6) 
,™ w w fl0\~ 1/4 / 7 \V2/10 14 GeV\ / m \3/2 

where g* is the number of relativistic degrees of freedom at temperature Trh and Mpi = 
1/V8ttG ~ 2.4 x 10 18 GeV is the reduced Planck mass. The reheating temperature must 
exceed about 5 MeV to preserve the predictions of nucleosynthesis [58l [59]. With m = 
10 3 GeV and 7 = 1, this puts an upper bound on v < 10 16 GeV, while for m = 200 GeV and 
7 = 0.1, the upper bound is strengthened to v < 10 14 GeV. We will mostly focus on values 
of the VEV less than v < 10 14 GeV for the rest of the paper. 

If flat-direction strings are to form, the corresponding U(l) gauge symmetry must be 
restored at or near the end of primordial inflation. Thus, if flat-direction strings are present 
in the universe today, they were most likely formed after a period of thermal inflation. The 
initial densities and properties of the strings depend on the details of the phase transition 
ending thermal inflation, when the flat-direction field overcomes the thermal barrier and 
starts to roll down to the true minimum. The nature of this transition has been studied 
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in Refs. [571 ITS"] . These authors find the tunnelling rate through the thermal barrier to 
be negligibly small until T ~ m. Below this temperature the tunnelling suppression is 
not parametrically large, and bubbles nucleate rapidly. Of particular importance to string 
formation is the radius of the bubbles of true vacuum when they coalesce, £. The initial 
size and separation between string segments are approximately equal to £. Since the phase 
transition proceeds quickly once the temperature falls below m, we expect £ to be within a 
few orders of magnitude of m -1 . 

The mechanism for string formation in the (a, b) model of flat-direction strings can be 
most easily understood in terms of flux-trapping. The winding number of a cosmic string is 
directly proportional to the net magnetic flux it carries in its core. In the broken phase, the 
magnetic flux is shielded. As a result, random fluctuations of the gauge field in the unbroken 
phase can be trapped between bubbles of broken phase. The scalar fields surrounding tubes 
of trapped flux then orient themselves to form a cosmic string with the appropriate flux 
quantum number. If £ is the typical bubble size at coalescence, the mean winding number 
of the strings formed in this way is [601 EI] 

N ~ y/tT f , (4.7) 

where Tf is the temperature at formation. Since the phase transition proceeds quickly once T 
falls below m, we expect that £ will not be too much larger than Tj ~ m -1 [15]. Therefore 
only the lowest winding modes will be significantly populated at the beginning. Let us also 
point out that the net magnetic flux of the configuration of Eq. (12.191) is N, independent of 
a and b. 

4.2 String Network Evolution 

Once cosmic strings are formed, their density evolves under the influence of the spacetime 
expansion, as well as the processes of reconnection and zippering. String reconnection is 
particularly important because it allows the string network to form loops and thereby transfer 
its energy into radiation. In the case of ordinary (abelian Higgs model) cosmic strings, the 
processes of string growth and loop production are found to balance each other, leading to 
a scaling solution. Flat-direction strings can also interact by zippering. This permits the 
formation of higher winding modes starting from an initial population consisting only of the 
lowest few modes. 

Cosmic string evolution has been studied extensively through numerical simulations [21 
[U [51 [6] . However, there has been no attempt that we know of to simulate a multi-tension 
string network including string zippering. In the absence of such simulations, we turn 
to analytic models of string evolution for guidance. A number of simple models have 
been constructed, and they give a good reproduction of the behavior of the long (horizon- 
length) string structure seen in simulations of the abelian Higgs model. To investigate the 
evolution of long flat-direction strings, we will make use of the model of Tye, Wyman, and 
Wasserman (TWW) [32], which generalizes the formulation of Ref. [7J. The TWW model 
was constructed to study the behavior of long superstring cosmic strings, which also exhibit 
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stable higher-winding modes and zippering, but with different rules for the outcome of string 
zippering. 

In the TWW model, long cosmic strings are characterized by a mean velocity u, a typical 
correlation length along the strings L, and a mean string number density n a , where a labels 
the winding number of the string (i.e. N = a). The number density of the string species a 
is defined through its relation to the energy density according to 

P° = n 2 ' ( 4 - 8 ) 



where \i a is the tension of the species. All string species are assumed to be described by 
the same v and L. This is a reasonable simplification for two reasons. First, the tension 
of different strings is a slowly varying function of the winding number, so in the absence of 
interactions with other string species, each string type should evolve in much the same way. 
Second, higher winding modes are mainly formed by the zippering of lower winding modes, 
and thus the speed and the fluctuation size of different string varieties should be roughly 
similar. 

The evolution equations for v and L in the TWW model are taken to be 

^ = HL(1 + v 2 ) + Cl u, (4.9) 



dv 
~dt 



'l-u 



2\ 



°j--v(2H) . (1.10) 



These equations are based on the model of Ref . [7] , where they are derived from the averaged 
equations of motion for a string evolving in an expanding Friedmann-Robertson- Walker 
background spacetimej§ 

The TWW model generalizes Ref. [7] by adding an independent density variable n a for 
each species. The value of n a is taken to evolve according to a Boltzmann-like equation 

h a = -2Hn a -- L j--P*n 2 a vL (4.11) 



+ FvL 



b,c 



^P abc n b n c (1 + 5 bc ) - P bca n c n a (1 + 5 ac ) 



Here, P a is proportional to the probability of self-reconnection for a string of variety a, 
P abc is the interaction probability for the process b + c — > a, and F is an overall non-self- 
interaction factor. Once the time dependence of H is specified, Eqs. (14.91 14.101 14. lip form a 
closed system describing the evolution of the long string component of a multi-tension string 
network. 

The values of the constants appearing in Eqs. (14.91 14.10114. IT]) can be fixed by comparing 
the scaling solution for a single (non-interacting) string to values obtained in string simula- 
tions. Ref. [32J reports that such an agreement is obtained with c\ = 0.21, C2 = 0.18, and 



6 Ref. [7] also considers frictional forces acting on cosmic strings. As in the TWW model, we do not 
include frictional effects in our analysis. We have checked that they are negligible for v > \/m Mpi, which is 
expected for the flat-directions strings under consideration. 
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Pi = 0.28. We use the same values for c\ and c% y which are related to the efficiency of loop 
formation and the amount of small-scale structure on the strings, respectively. For P a and 
F, we set them to P a = F = 0.28/2 = 0.14. Since P a is proportional to the probability of 
reconnection, this accounts for our assumption that a pair of strings is just as likely to zipper 
as to reconnect when both outcomes are kinematically allowed. We also set the coefficients 
P a bc to 

p = f 1; a=\b±c\, v < z/ thresh , 
abc \ 0; otherwise. 

These values are in accord with our assumptions about zippering. Motivated by the results 
of Section [31 we set the velocity threshold for zippering to thresh = 0.85 in our numerical 
analysis. 

To evaluate Eqs. (14.91 14.101 14. lip describing the evolution of the string network, we must 
also specify the evolution of the Hubble parameter H appearing in these equations. We 
do this by solving for the scale factor a(t) using Eq. ( 14. 4ft . After thermal inflation, the two 
dominant sources of energy density are p^,, from the oscillations of the light scalar field, and p r 
for radiation. We begin the evolution at the time U = 10 Mp\/mv, as would be expected after 
thermal inflation. The initial radiation density is taken to be p r (U) = m 4 , while the initial 
scalar field energy density is set to p^iti) = m 2 v 2 . After time U, p^ and p r evolve according 
to Eq. ( I4.2p . Since we are interested in running the string evolution equations all the way to 
the present time, we also add a very small matter density at the end of thermal inflation, 
at tRH = r^ 1 . The initial matter density is chosen such that it becomes the dominant form 
of energy at the approximate equality time t eq = 3 x 10 36 GeV" 1 . For reference, the present 
time is about t ~ 6.6 x 10 41 GeV -1 . With m = 10 3 GeV, v = 10 13 GeV, and 7 = 0.1, the 
initial matter density is PmiiRu) = (8.0 x 10~ 3 GeV) 4 . At later times, this dilutes according 
to dp m jdt = —3Hp m . Throughout the evolution of H, we self-consistently assume that the 
energy density due to the string network plays a negligible role. 

We appeal to our expectations from thermal inflation to set the initial values of the 
variables v, L, and n a . The symmetry breaking phase transition after thermal inflation 
occurs quickly once the temperature falls below T = m. The mean bubble radius £ when 
they coalesce should be therefore not much larger than the nucleation radius, which is close 
to m -1 [15]. Thus, we set = 5m" 1 and ni(tj) = l/(5m -1 ) 2 as reasonable starting values. 
The initial densities of the higher winding modes, a > 1, are set to zero. We also choose 
v(ti) = 0.9. While there is considerable arbitrariness in these choices of initial conditions, 
we find that our results at late times are largely independent of them. 

In Figs. [10] and [UJ we show the numerical solutions of the string network equations for 
the model parameter values m = 10 3 GeV, v = 10 13 GeV, and 7 = 0.1. For comparison 
with Section [21 this choice corresponds to a value of A = g 2 m 2 /v 2 ~ 10 -20 . Fig. [TUI depicts 
the evolution of the densities of the five lowest winding modes in terms of the quantities 

^ = ^ a = -^P=, (4.13) 



Pa PcVT 



where Q a is the ratio of the energy density of string species a relative to the critical density 
p c = 3 H 2 /87rG, and p a is the tension of string species a. Normalizing by the tension makes 
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Figure 10: Evolution of cosmic string densities after thermal inflation with v = 10 13 GeV, 
m = 10 3 GeV, and 7 = 0.1. We have also set N max = 50 in generating this plot. 

fl a proportional to n a times a quantity that is independent of the winding number. In Fig. [11] 
we show the evolution of the universal length scale L and universal string velocity v. 

Figs. [TU1 and [TT1 show that (within the TWW model) the string energy densities approach 
a scaling solution at late times as evidenced by HL, v, and Cl a all flowing to constant values. 
The scaling length, velocity, and densities are largely independent of the initial state of the 
string network. At late times, the string densities make up a nearly fixed fraction of the total 
energy density of the universe. We also find that the early era of oscillation dominance during 
reheating does not alter the final string densities in an appreciable way. These features are 
very similar to what is found in simulations of ordinary (abelian Higgs) string networks with 
only a single string species [5] [6] . 

The interesting new feature in the evolution of flat-direction cosmic strings is that nearly 
all string species flow towards very similar scaling values. This is the result of string zippering, 
which allows the formation of higher winding modes from lower ones. Note that the formation 
of these higher modes does not begin immediately. With the initial values specified above, 
the initial string length scale L is much smaller than its scaling value, which is close to the 
horizon scale. This has the effect of rapidly driving the string speed to its maximal value, 
v — > 1, at the outset, as can be seen in Fig. [TTJ, which effectively shuts off string zippering. 
Once L and v settle down to near their scaling values, zippering begins and the higher 
winding-mode densities quickly flow towards their scaling values. This scaling behavior is 
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Figure 11: Evolution of cosmic string speed and length scales in the aftermath of thermal 
inflation with v = 10 13 GeV, m = 10 3 GeV, and 7 = 0.1. We have also set N max = 50 in 
generating this plot. 

quite robust. Changing the values of F and z/ thresh does not alter the qualitative string 
densities provided ^thresh is larger than the mean string velocity in the scaling regime. 

The fact that many string species flow towards equal scaling values complicates the 
numerical analysis, since numerical limitations allow us to include only a finite number of 
winding modes up to an unphysical maximal value N max . In making Figs. [10] and [11] we have 
set N max = 50. We also find that the final, nearly universal scaling density of the strings 
depends on the artificial value of N max . This feature is illustrated in Fig. [12] To a good 
approximation, the near-universal string scaling density goes like 

&a oc (4.14) 

1 \ rr 



' max 



as illustrated by the dotted line in Fig. [12] Evidently, the string energy density gets spread 
out among the many string types. There is also the question of how to handle the zippering 
of strings whose winding numbers sum to greater than N max . In principle, these strings can 
zipper into modes with iV > N max which are not included in the simulation. In Figs. [10] 
and [11] and in the analyses to follow, we simply disallow all such zippering processes. This 
leads to slight increase in the scaling density of modes with N > N max /2. However, we have 
also studied other prescriptions for handling these zippering events, and for the examples we 
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Figure 12: Dependence of the scaling-regime string density on the total number of string 
species included in the simulation, N max . The dotted line shows a fit to O x oc 1/N max . 

looked at, we find qualitatively similar results for the modes with iV -C N max . 

The dependence of the scaling densities on N max is clearly unphysical. We would like 
to take N max —>■ oo, but this has its own problems. Since the energy density at large N 
goes like In a (from the logarithmic dependence of the tension on the winding number), if all 
string species flow towards a universal scaling density proportional to N max the total network 
energy density goes like 

Ptot oc — ^ lna ~ \nN max . (4.15) 

max , 
a=l 

This diverges logarithmically as N max — > oo. In practice, however, this divergence is not 
realized. The initial string spectrum consists almost entirely of the lowest modes, the density 
of higher modes is built up from the lower modes by zippering, and these higher modes take 
longer to reach their scaling values. At any given time, only a finite number of strings have 
developed their scaling density.0 Let us define N eq (t) as the highest mode that has reached 
scaling by time t. Modes with N > N eq (t) all have densities well below their equilibrium 
scaling values. Thus, at time t, we effectively have N max = N eq (t), and the total energy 
contained in the string network goes like In N eq (t). 

In Fig. [T3] we show the time evolution of N eq (t) for several values of N max . All other 
parameters are the same as in Figs. [TUJ The curves for different values of N max match up 



7 In this sense, our use of the term scaling for flat strings is somewhat more general than its meaning for 
ordinary cosmic strings because the string densities are not completely static, but very slowly varying. 
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for iV < N max /3, but start to deviate from each other as the winding number N approaches 
Nmax- Focusing on the apparently universal portion of these curves, the rate of increase of 
N eq (t) with time goes like t 0,22 . If we can extrapolate this dependence to much larger winding 
numbers, the value of N eq at the present time to will be 



where we have used tj ~ 10 M P i/mv ~ 2.4 x 10 3 GeV _1 and t ^ 6.6 x 10 41 GeV -1 . This is 
a very large number, but it is not so large so as to be problematic. Recall that the string 
tension, given in Eq. (I2.49p . increases logarithmically with the winding number. The tension 
of a string with N = 10 8 is merely 



for m = 10 3 GeV and v = 10 13 GeV, corresponding to A ~ 10~ 20 in Eq. (12.491) . Moreover, 
the total string energy density in the network is less than about \iaN eq (t ) < 20 times the 
energy density of a network containing a single type of string with the same tension as the 
lowest mode. These values for the maximal tension and the total string density are not 
much larger than for an ordinary cosmic string, and they present no obvious cosmological 
difficulties. 

In our analysis of flat-direction string network evolution described above we have used a 
very simple analytic model of string network evolution; we have made specific assumptions 
about the details of the string interactions; and, we have made extrapolations into regions 
well beyond what we are able to probe analytically and numerically. Nevertheless, a simple 
picture for the evolution of a flat-direction string network emerges from our results, and is 
likely to be genuine, even if some of the underlying assumptions are not necessarily rigorous 
and the model used to study the network evolution is overly simple. In this picture, a very 
large number of string species reach similar scaling densities by the present time. The total 
energy of the network is within an order of magnitude or two of the energy density that 
a single abelian Higgs string species would have for the same value of the string tension. 
However, instead of being concentrated within a single species, the string energy density is 
nearly uniformly distributed among all the string species that have attained scaling. Thus, 
the flat-direction string network consists of a near continuum of string species, but with 
global properties that closely resemble those of a single species network. 

Our argument for this picture is based on the very slow dependence of the flat-direction 
string tensions on the winding number. On account of this slow variation, the macroscopic 
properties of the many species that have attained scaling are very similar to each other. For 
example, Fig. [2] shows that the tension of a mode with N = 100 is only about 1.4 times that 
of the N = 1 mode for A ~ 10~ 20 . From this feature, as long as the zippering is reasonably 
efficient and the lowest mode is able to attain a scaling value for its density, we expect the 
densities of the string modes to be very similar to one another up to large values of the 
winding number N ^> 1. 

One curious aspect of this picture is that the total energy density in the network corre- 
sponds to less than a few hundred individual strings of horizon length. It is therefore curious 




(4.16) 



[L N < 3yUl, 



(4.17) 
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Figure 13: The number of string species that have reached scaling, N eq , as a function of 
time in the aftermath of thermal inflation for different values of N max , with the parameter 
values v = 10 13 GeV and m = 10 3 GeV. The dotted line indicates an approximate fit to 
N eq (t) oc t a22 in the region where the curves appear to be universal. 

that the TWW model applied to flat-direction strings predicts that there are many more 
string species than this in the scaling regime at the present time, each with a characteristic 
length scale of horizon size. We suggest that the scaling densities predicted by the TWW 
model for flat-direction strings should be interpreted as time-averaged values. At any given 
epoch in the scaling regime, there exist many fewer long strings than N eq (t). However, these 
strings are continually zippering into other string species, and averaged over time, many 
more string species are populated (with a lower density per string) than are present at any 
one time. It is also possible that this issue of discreteness leads to a value of N eq (t) that is 
smaller than what is predicted by the TWW model. 

A definite confirmation of this picture of flat-direction string evolution would appear to 
require a full numerical simulation of the network (as well as lattice simulations to determine 
the zippering probabilities). This task is complicated by the need to include many different 
string species in the simulation, and is beyond the scope of this paper. We have, however, 
examined the effect of changing some of our assumptions about string zippering encoded in 
the coefficients P a bc, defined in Eq. (14.121) . For example, we find that reducing the probability 
for zippering into a = (b + c) relative to a = \b — c| does not significantly alter the final 
scaling densities. We have also looked into modifying the interaction terms in Eq. (14. lip . 



34 



as suggested in Ref. [35], and we again find the same qualitative picture of string network 
evolution. These results suggest that the picture of flat string evolution presented here is 
robust. 



Before moving on, let us briefly compare our near- continuum picture of flat-direction 
cosmic strings to the cosmological picture of (p, q) surperstring cosmic strings derived in 
Ref. [52] . These quasi-fundamental strings can be labelled by pairs of integers (p, q) with 
p > 0. A string state is stable only if p and q are relatively prime. States with (p, q) not 
relatively prime can be formed but are only marginally stable. They are expected to decay 
into lower, stable modes after they are created. In the analysis of Ref. [32J , this additional 
dissipative channel led to a rapid decrease in the relative population of higher-tension modes. 
That superstring cosmic strings do not form a near-continuum scaling network is also not 
surprising given that the tensions of these strings increase fairly rapidly with the mode 
numbers 132 1, 



where g s is the superstring couplingjfl Hence, even though flat-direction cosmic strings and 
(p, q) strings can both form stable winding modes through zippering, these two varieties of 
cosmic strings interact and evolve in significantly different ways. 

5 String Signatures 

If cosmic strings are present in the early universe they can give rise to a number of observable 
signatures. No evidence for cosmic strings has been found in the temperature power spectrum 
of the cosmic microwave background or in large-scale sky surveys. This implies the constraint 
C*H < 3 x 10~ 7 , fairly independently of the underlying string model [63J, EU [65j [66j 157] . 
Beyond these limits, the most promising signatures for ordinary (abelian Higgs) cosmic 
strings are gravitational lensing and gravitational radiation [2113]. We find that these signals 
can be modified for flat-direction cosmic strings. Flat-direction cosmic strings are also more 
likely to radiate into their constituent particles than ordinary cosmic strings, leading to new 
classes of potential signatures. By combining observations of several different phenomena, it 
may be possible to distinguish flat-direction cosmic strings from ordinary cosmic strings as 
well as (p, q) cosmic superstrings. 

5.1 Gravitational and Particle Radiation from Loops 

Cosmic strings emit gravitational radiation primarily through the oscillations of string loops. 
For both ordinary and flat-direction cosmic strings, a single loop is expected to emit gravi- 
tational radiation with power 



This formula applies in ten-dimensional flat space. It may receive corrections in other backgrounds [62j . 
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where T = 10—100 is a dimensionless constant whose precise value depends on how the loop 
is oscillating [691 CHS EJ E21 This rate is independent of the length of the loop, I. The 
radiation frequencies do depend on I and are 



n 



2n 



n= 1,2,3,... 



(5.2) 



with the relative power going into mode n decreasing at least as quickly as n ' for simple 
string loop solutions [701 17T1 1721 173] . 

To compute the gravitational wave background from a cosmic string network, one must 
convolute the power emitted by individual loops with the loop density distribution. Unfortu- 
nately, even for ordinary cosmic strings, the loop density distribution is not fully understood. 
The main uncertainty is the size of loops when they are formed. It is standard to parametrize 
the typical initial loop length according to 



where t is the time of loop formation, and estimates for a range between the string width [8], 
to (r Gfi) x with x > 1 [ZU [751 [76], all the way up to a = 0.1 [77] . We will consider different 
values of a below. 

Cosmic string loops can also radiate directly into particles [78j [79j [80]. This can arise 
both through the direct emission of particles from smooth strings [78j [791 [80] , as well as from 
cusp annihilation [81]. For the string loops present in the early universe, cusp annihilation 
is usually the more important source of particle emission [81] A cusp is a point on a string 
that reaches the speed of light at some instant during its (Nambu-Goto) evolution. Cusps 
are a generic feature of many simple solutions for the motion of a string loop, where they are 
found to occur about once per oscillation period J70j[3]. In the region near the cusp, the string 
segments fold back upon themselves such that the separation between the adjacent segments 
becomes smaller than the string width. This allows these string segments to annihilate each 
other. Cusps should not be confused with string kinks, which are points on a string where the 
tangent vector changes substantially over a very short distance, on the order of the string 
width [73]. Unlike at a cusp, there need not be any significant annihilation of the string 
segments in the vicinity of a kink, and kinks can persist for many loop oscillations [83J. 
Kinks can be created from string reconnection and zippering. 

The effective length of the overlap region between the adjacent string segments near a 
string cusp on a loop of length £ is about 



where w is the string width [84J. The overlapping string segments near the cusp are 
expected to annihilate, transferring most of the string energy within the overlap region 

9 This conclusion can change if there exist light (superstring) moduli fields with masses much smaller 
than where w is the width of the string [79l [80] . For flat direction strings, both the string width and 
the typical moduli mass are set by the scale of supersymmetry breaking m. As a result, the rate of moduli 
emission by flat-direction strings is suppressed, and the corresponding bounds |82j are not relevant. 



ii = a t, 



(5.3) 




(5.4) 
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to the constituent particles making up the string. The total average power released into 
particles through this process by a single string loop is [M] 



Pcusp — f^^c ( „ ) • (5-5) 



Here, cjl is the cusp rate, where l/£ corresponds to the period of a loop oscillation, and c 
is the probability per period for a cusp to occur. We expect c ~ 1, although it has been 
argued that the presence of kinks on strings could push it to smaller values [73J. This is yet 
another uncertainty associated with the structure of cosmic strings on small scales. 

For a given tension, flat-direction strings are much wider than ordinary cosmic strings; 
w ~ m~ l v~ l compared to w ~ v ~ 1 ~ ji~ l l 2 . The amount of string annihilated in a cusp 
is therefore greatly enhanced. The total particle radiation power from cusp annihilation by 
a flat-direction string loop is 

Pcusp — / (5-6) 



Relative to the gravitational radiation power, Eq. ( 15. II) . we see that cusp annihilation 
dominates for sufficiently small loop sizes. The loop size at which the two powers become 
equal is 



u * m {TGt)- (57) 

Recall that the loop size at formation is £j = a t. For l± < £ = the loops will decay primarily 
through particle emission, and not gravitational radiation. On the other hand, when ^ 3> £=, 
most of the loop energy will go into gravity waves, except for a small burst of particles towards 
the end of the loop's existence. 

Thus, the emission of particles by flat-direction cosmic strings through cusp annihilation 
is greatly enhanced relative to ordinary cosmic strings. If cusp annihilation dominates over 
gravitational radiation, many of the gravitational radiation signals will be suppressed com- 
pared to ordinary cosmic strings. In order to compare the relative signals from gravitational 
radiation and particle emission, it is helpful to concentrate on three particular epochs in the 
early universe: the reheating time the time at which £^ — t = — at = ; and the earliest 
time tf at which a given gravitational wave frequency mode / can form. 

We found in Section H] that reheating after thermal inflation occurs when t ~ £r# := T^ 1 , 
where = •ym 3 /v 2 is the decay rate of the light flat-direction scalar field. This yields 

/0.1\ / V \ 2 flO 3 GeVV/^20 r, ,r IN /r on 

Recall that if the process of reheating after thermal inflation is to avoid disturbing the 
predictions of nucleosynthesis, we must have v < 10 16 GeV for m = 1000 GeV and 7 = 1, 
and v < 10 14 GeV for m = 200 GeV and 7 = 0.1. 
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The second moment of interest, the time after which newly-formed loops lose most of 
their energy in the form of gravity waves, occurs when £j = l = = at = . This corresponds to 
the time 



~ a 1 c 2 



50\ 2 / 2x lO^ y AO 3 GeV \ 



(10 15 GeV- 1 ). 



We have expressed t= in terms of G[i rather than the VEV v because it is this dimensionless 
combination that appears frequently in the estimates below. An approximate conversion 
between G\i and v is (see Fig. CP 

Gt^( ° S (2X10- 11 ). (5.10) 



10 14 GeV 

Given the upper bound on v from reheating after thermal inflation, we will mostly focus on 
v < 10 14 GeV. 

The third time of interest is tf, the earliest moment at which a given gravitational wave 
frequency as low as / can be emitted. Recall that loops formed at time have the initial 
size £{ti) = ati, and subsequently shrink and radiate into frequencies / > 2/E. For a mode 
observed at the present time with frequency / = f{to) emitted at time t, the initial frequency 

was 

a[t) 

Combining these facts, the earliest time tf at which a mode with present frequency / could 
have been emitted is 

af a(t ) 



( _/__„„ \ 3 

a 



-3 fw-THz\ (a5 x 1Q 10 GeV -l) tf > te 



a~ 2 (^j 1 ^) (2.5 x 10 19 GeV" 1 ) t RH < t f < t eq 

a" 3 (^) 3 ('-^f 1 ) 112 (1-2 x 10 22 GeV" 1 ) tf < t RH 



Both t= and tf depend on the parameter a that characterizes the typical size of a string 
loop when it is formed, £i = ati. The dynamics of loop formation are not completely 
understood, and as a result, estimates for a vary widely. Some recent simulations find that a 
significant portion of the loops formed are quite large, with a ~ 0.001 [851 EE] or a — 0.1 [77] . 
Other simulations find that the typical initial loop size approaches their resolution limits [8]. 
In this case, it is thought that gravitational radiation will smooth out very small fluctuations, 
and impose a lower limit on a (3j. The scale over which this smoothing occurs is also under 
ongoing investigation. Early estimates suggested a = T G[i [3] , but more recent analyses 
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have found even smaller values of a. In Ref. [73] the authors obtain a = (TGfi) x with 
X = 1-5 during the radiation era and x = 2.5 during matter dominance. The authors of 
Ref. [75] find a ~ 0.6 T {G^i) x with % = 1.2 in the radiation era and x = 1-5 in the matter 
era. Furthermore, in Ref. [76] it is suggested that the simulation results of Ref. [77] should 
be interpreted as predicting a network with 10-20% of the loop energy density in the form of 
large loops with a ~ 0.1 and the remainder in the form of very small loops with a ~ T(Gfj,) x 
with x > 1- O n account of the rapidly evolving state of the field, we will consider both large 
and small values of a below. 

Fig. [TJ] shows tun, t=, and tf in GeV -1 units as functions of the VEV v for large loops 
with a = 0.1. The model parameters were set to m = 10 3 GeV, c = 1, 7 = 0.1, T = 50. 
In this plot we also indicate the present time to — 6.6 x 10 41 GeV and the matter-radiation 
equality time t eq ~ 3.5 x 10 36 GeV with dotted lines. The value of tf is shown for two values 
of the frequency, / = 10~ 7 Hz and / = 10 2 Hz. These values span most of the range relevant 
for gravitational wave searches. For these large loops, t= is always much less than t eq , and 
all loops formed after t= will decay predominantly into gravitational radiation. At the lower 
frequency / = 10~ 7 Hz, tf lies below t eq but above tim, and is never much less than t=. This 
suggests that the gravitational wave signal at this frequency will not be attenuated much by 
the enhanced rate of particle emission by the loops. On the other hand, the value of tf for 
/ = 10 2 Hz lies well below both t^H and t = , indicating that the high-frequency gravitational 
wave signal will be reduced. 

The values of tnH, t=, and tf for very small loops, a = 0.6 T (G/x) 1,5 , are shown in Fig. [15] 
as a function of the VEV v. This value of a corresponds to the estimate of Ref. [75] for loops 
emitted in the matter era. Even smaller values of a are suggested in Ref. [74J. As before, 
the other model parameters were taken to be m = 10 3 GeV, c = 1, 7 = 0.1, and T = 50, 
and the dotted lines denote the present time to — 6.6 x 10 41 GeV and the matter-radiation 
equality time t eq ~ 3.5 x 10 36 GeV. This figure indicates that the prospects for gravitational 
radiation from small flat-direction string loops are much less promising than for large loops. 
Indeed, t= is larger than the present time to for v < 2 x 10 12 GeV. Such small loops will 
decay almost entirely to particles instead of gravitational radiation. Even when t = is less 
than to, the curves for tf show that the gravitational wave signal is very suppressed relative 
to the signal from large loops. At low frequencies / ~ 10~ 7 Hz there is no signal at all since 
tf exceeds the present time to', the loops are simply too small to radiate into this frequency 
range. Even for frequencies near / = 10 2 Hz, there will be a gravitational wave signal only 
for v > 3 x 10 12 GeV. Despite the reduction in the gravitational wave signal, small loops 
may be observable through their copious emission of particles. 

In summary, we find that the large width of flat-direction cosmic strings greatly enhances 
the rate at which they decay into their constituent particles through cusp annihilation. With 
this enhancement, our preliminary analysis indicates that string loops that are initially large 
(a ~ 0.1) decay predominantly into gravitational waves, while very small loops (a rG/i) 
decay primarily into particles. The typical size of string loops when they are formed is an 
unresolved problem, and well-motivated arguments in favor of large loops, very small loops, 
or possibly both at once, can be found in the literature. In the face of this uncertainty, we will 
focus on two particular choices of the loop size parameter a to estimate the observational 
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Figure 14: Dependence of the times of interest tjuj, t = , and tf on the VEV v for a = 0.1. The 
black dotted lines indicate the present time to — 6.6 x 10 41 GeV -1 and the matter-radiation 
equality time t eq ~ 3.5 x 10 36 GeV -1 . 

signatures from flat-direction cosmic strings. To compute the gravitational wave signals 
we will set a = 0.1 for all loops, as suggested in Ref. [77]. Our results can be rescaled 
appropriately when only a fraction of the loops are large. To estimate the signals from 
particle emission due to cusp annihilation, we will instead assume that a is sufficiently small 
that all loops decay mostly into particles. This is plausible for flat-direction strings for which 
the rate of particle emission by cusp annihilation is enhanced. Again, it is straightforward 
to modify our results to accommodate larger values of a. 

Finally, let us also mention that the picture of loop formation by flat-direction strings 
might be different from that of ordinary cosmic strings. For example, the enhanced rate 
of particle emission by cusp annihilation could potentially smooth out small fluctuations 
over scales larger than the (naive) gravitational radiation scale TGnt. String loops can also 
remain bound to the parent string, as illustrated in Fig. [HJ This could modify the distribution 
of initial string loop sizes. The rate of cusp formation on these bound loops may also be 
different from that on free loops. 
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Figure 15: Dependence of the times of interest tjm, t=, and tf on the VEV v for the 
representative small initial loop size parameter a = 0.6 T (Gfj,) 1 ' 5 . The black dotted lines 
indicate the present time to ~ 6.6 x 10 41 GeV -1 and the matter-radiation equality time 
t pa ~ 3.5 x 10 36 GeV" 1 . 



5.2 Gravitational Wave Signatures 



Cosmic strings can give rise to two types of gravitational wave signals. The combination 
of many string loop decays produces a smooth stochastic background of gravitational radia- 
tion [69J. On top of this background, individual cusps can produce intense bursts of gravity 
waves [87]. Gravitational wave detectors are sensitive to both types of signals. For the 
string tensions of interest, G/i < 1CT 10 , the stochastic background is the more promising 
one [88j [89j EQ] and we will focus on it. To estimate this gravitational wave background due 
to flat-direction cosmic strings we will assume that all string loops are large when they are 
formed, with a ~ 0.1 [77]. If only a fraction of the loops produced are large, as advocated 
in Ref. [76], our results can be rescaled by this fraction. 

We compute the gravitational radiation density due to cosmic string decays following 
Ref. [72]. Consider radiation in the frequency range (/, / + df) observed today that was 
emitted at time t. Keeping track of only the lowest mode0 this radiation was emitted by 

10 In Ref. [91] this was shown to be a good approximation for computing the stochastic background. 
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loops of size {£ — d£, £), where 
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These relations follow from the loop evolution equation £(t) = otU — TGp(t — ti), valid for 
t > U and £{t) > £=. 

The rate at which loops are formed during the string scaling regime can be estimated using 
the results of numerical simulations or from simple analytic models like the one presented in 
Section HI These predict a net energy flux into loops of 
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with ( ~ 10, and being the scaling energy density of long strings. This result can be 
obtained by summing Eq. (14.111) over all string species that have equilibrated. It follows that 
the rate per unit volume that loops of initial size a t are formed is 
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Applying this result to loops formed in the time range (ti — dti,ti), the number density 
of loops radiating into the frequency range of interest at time t is 
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The redshift factor in this expression accounts for the dilution of the loops as they evolve 
from ti to t. Given that each loop radiates gravity waves with a power TGp 2 , we can combine 
everything and sum over t to find the signal. The total gravitational wave density at the 
present frequency / is 
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Figure 16: Gravitational wave density for flat-direction cosmic strings as a function of 
frequency for four different values of the VEV v. The solid lines include the cutoff I > £= 
due to cusp annihilation. The dashed lines show what the gravitational wave density would 
be without this cutoff. 

Here, p c is the critical density, and ti and £ are functions of t defined by Eqs. (15.131) and (I5.14p . 
The integration limits range between if := max(tf, 10 5 GeV -1 ) and to, where tf is given in 
Eq. (15.121) 1"] Noting that a oc t 2 ^ 3 during the matter era (t < tuji and t > t eq ) and a oc t 1 ^ 2 
during the radiation era < t < t eq ), this equation can be integrated straightforwardly. 
Relative to the treatment of Ref . [72] , we have included a cutoff of £ > £= = a t = through a 
step function. This accounts for the loops only being able to radiate efficiently into gravity 
waves if their length is greater than £ = . It is this cutoff, along with the additional redshifting 
that occurs during reheating after thermal inflation, that suppresses the gravitational wave 
signal from flat-direction strings compared to ordinary strings. 

In Fig. [11] we show the stochastic gravitational wave signal from initially large cosmic 
string loops as a function of frequency. We have used the parameter values a = 0.1, T = 50, 
m = 10 3 GeV, and 7 = 0.1 in making this plot. The solid lines show the gravitational wave 
density from flat-direction cosmic strings, including the cutoff £ > £ = and the additional 

^^Normally the lower limit would simply be tf, but in the present case the flat-direction string network 
only reaches scaling at t ~ 10 5 GeV -1 . Numerically, we find that this additional cutoff has no visible effect 
because the gravity waves emitted shortly after the end of thermal inflation are diluted away during the 
subsequent reheating. Gravity waves from the phase transition [92 will also be diluted by thermal inflation. 
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redshifting during reheating after thermal inflation. The dashed lines indicate what the 
signal would be for ordinary cosmic strings, without the cutoff i > £= or reheating effects. 
At lower frequencies the relevant loops are formed later on, at times greater than tjnj an d 
t = , and there is no change to the signal. At higher frequencies, the cutoff on the loop size 
and the additional dilution during reheating both suppress the gravitational wave signals. 
As can be seen in Fig. [T3J the cutoff i > £= = a t= is more important for lower values of v 
(and G/j,), while the reheating dilution is more significant at larger values of v since tmj is 
larger. This is why the shape of the high frequency cutoff changes as we increase v. 

The attenuation of high frequency gravitational wave signals is relevant to LIGO and 
Advanced LIGO, which can potentially probe down to Qcwif) — 10~ 9 at frequencies around 
/ = 10 2 Hz [93J. Fig. [16] indicates that LIGO is not expected to be able to find evidence for 
flat-direction cosmic strings. On the other hand, the prospects for discovery at LISA and 
from measurements of pulsar timing are quite promising. The LISA probe is expected to 
cover portions of the range 10~ 4 Hz < / < 1CT 2 Hz down to ^Gw(f) — 10~ n [91]. Since 
the gravitational wave signal from large flat-direction string loops is mostly unmodified in 
this frequency range, LISA will be able to probe a sizeable portion of the model parameter 
space. Limits from pulsar timing are currently ^Gw{f) ^ 3 x 1CT 8 in the frequency range 
1CT 7 — 10~ 8 Hz [95], which is again low enough that the gravitational wave signal from 
flat-direction cosmic strings is unsuppressed. From this bound we obtain the constraint 
v < 10 14 GeV. It is expected that this limit will be improved to Slowif) ~ 1CT 10 by 
upcoming experiments [95]. Note that flat-direction cosmic strings offer the interesting 
possibility that LISA and pulsar timing experiments could detect a stochastic gravitational 
wave background with Vt GW (f) > 10~ 9 , while (Advanced) LIGO sees nothing even though it 
is sensitive to signals at this level. This would be a suggestive hint for flat-direction cosmic 
strings. 

5.3 Particle Emission Signatures: Dark Matter 

Having studied the gravitational wave signatures of string loops that are large when they 
are formed, let us now consider the possibility that the typical initial loop size is very small, 
aCT G/i, as suggested in Refs. [7H [75]. If a is sufficiently small, nearly all the energy of a 
loop is released as particle excitations of the fields making up the string. This is plausible 
for flat-direction cosmic strings due to their enhanced rate of particle emission by cusp 
annihilation relative to ordinary cosmic strings. In the (a, b) models of flat-direction strings 
presented in Section [2j the fields making up the string consist of two chiral supermultiplets 
and one massless vector (gauge) supermultiplet. When the chiral supermultiplets develop 
VEVs, it is more convenient to describe the theory in terms of a heavy massive vector 
supermultiplet, with mass on the order of gv, as well as a light supermultiplet with mass on 
the order of m [44J. This light multiplet is light on account of the flatness of the potential. 
Cusp annihilation will produce both the heavy and the light states making up the string. 
These particles will subsequently decay, and can be a potential source of dark matter and 
high energy cosmic rays. We consider both of these possible signatures in turn, assuming 
that all string loops are very small and decay entirely into particles rather than gravitational 
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waves. 



To compute the dark matter density from decaying string loops we again make use of 
Eq. (I5.16P that specifies the rate at which the scaling string network transfers its energy into 
loops. Contributions to the dark matter density from loops produced before the network 
attains scaling are diluted away by the subsequent reheating process^! We assume that 
a fraction e% of the energy emitted by the cusp annihilations of loops eventually becomes 
dark matter (such neutralino or gravitino LSP)@ The total dark matter density at the 
present time from the strings is then 
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where to is the present time and tf is the time at which the DM particles freeze out 
of equilibrium. The factor of [a(t)/a ] 3 accounts for the additional dilution of the dark 
matter (or the constituent string fields) after they are produced. It is convenient to split the 
integration into three pieces: t eq < t, tun < t < t eq , and tf <t < £r#. These integrations 
are straightforward and yield 
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where 7 is the prefactor appearing in Eq. (14.51) . Numerically, the largest contribution comes 
from the third term, from the integration range tf Q < t < tRH. We have bounded the 
logarithm in this term from above in making this estimate. For reasonable values of the 
model parameters, the amount of dark matter produced by decaying loops is safely small, 
although smaller values of v are preferred. This differs from the much stronger constraints 
on regular cosmic strings that are able to decay into dark matter [SB], which is due to the 
dilution from reheating after thermal inflation. 

There is an additional contribution to the DM density from the out-of-equilibrium decays 
of the oscillating flat-direction fields during reheating. If a small fraction €2 of these decays 
ends up as dark matter, the present contribution to the DM energy density will be 
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12 We have verified this using the simulation of Section [4] 

13 The decay products from a cusp annihilation are typically boosted by factor of (i/w) 1 ^ 2 [53] , This 
can lead to values of e\ <^ 1 if most of the cusp energy goes to kinetic energy rather than non-relativistic 
dark matter. For small initial loop sizes, a < TGfi, this boost is weak enough that the dark matter decay 
products quickly reach kinetic equilibrium for T > Tru 1 when most of the dark matter is produced. 
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Thus, the branching fraction into DM particles must be very small. Note that e\ and 62 
can be very different from each other. The particles emitted from a cusp annihilation can 
include some of the heavier component fields making up the string. On the other hand, the 
decays of the flat-direction fields after thermal inflation involve only the light modes. The 
decays of these states into superpartners (such as a neutralino or heavier gravitino LSP) can 
therefore be highly suppressed or kinematically inaccessible, allowing for e<i <C ei- 



5.4 Particle Emission Signatures: Visible Matter 

In addition to dark matter, the decays of very small string loops can produce hadrons, 
leptons, and photons. This particle injection will be spread out over time as the scaling 
string network continually rids itself of excess energy by forming loops. Visible particles 
created by loop decays can imprint themselves upon the early universe in a number of ways. 
The energetic products from loop decays at temperatures below 5 MeV can disrupt the 
predictions of big-bang nucleosynthesis (BBN). At later times, energetic photons from loop 
decays can modify the blackbody spectrum of the cosmic microwave background (CMB). 
Some of the decay products from string loops can also be highly energetic, producing ultra- 
high-energy cosmic rays and contributing to the extragalactic diffuse gamma-ray background. 
We consider the possible signatures from cosmic strings from each of these effects. As for 
our dark matter estimates, we assume that all loops are so small that they decay entirely 
into particles. 

To estimate the effects of decaying string loops on BBN we make use of the results of 
Ref. [97J. In this work the authors used the successful predictions of BBN to place limits 
on the combination mxYx for a long-lived particle X of mass mx, relic density (per unit 
entropy) Yx, and lifetime tx, decaying at time t ~ tx- In contrast to a long-lived relic particle 
whose decays can be treated as being instantaneous, cosmic string loops are produced and 
decay continuously. These decays therefore have a cumulative effect on the light element 
abundances. To obtain a limit for decaying string loops we interpret the bounds from 
Ref. [97J as limits on the total energy injected within a comoving volume, mxYx = AE/S, 
where S is the total entropy within the volume a 3 . 

The total energy injected into the comoving volume a 3 by string loops that decay during 
the time interval (t a , t&) is 
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In writing this expression, we have implicitly assumed that t a and 4 both lie within the 
era of radiation dominance, as is relevant for BBN. The strongest limits on energy injection 
from BBN come from the relative fractions of deuterium and lithium-6 relative to hydrogen. 
Both of these are formed at times later than t > 100 s. Since the visible decay products 
from the loops thermalize quickly relative to the Hubble time, we set t a = 100 s to find 
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the bounds due to the deuterium and lithium-6 abundances [97]. Assuming a hadronic 
branching fraction of order unity, the total energy injection per unit entropy must be less 
than AE/S < KT 14 GeV0 This bound is satisfied provided v < 10 13 GeV. 

Late-time energy injection is also constrained by the the nearly perfect blackbody spec- 
trum of the CMB observed by COBE/FIRAS [98J. Photons produced by the decays of 
string loops that occur after the time tdc — 10 31 GeV -1 can distort this spectrum. Before 
tac, double Compton scattering (e + 7— > e + 7 + 7) efficiently thermalizes any additional 
photons that are created. The precise form of the spectral distortions created after t dC 
depends on the time at which the photons were injected. However, the net constraint from 
the non-observation of such distortions can be reduced to a constraint on the total photon 
energy created after tdc, Ap 7 /p 7 < 7 x 1CT 5 [99, 11001 1101] . The net photon injection from 
decaying string loops can be estimated using the rate of energy deposition by the network. If 
all the energy injected is in the form of photons (possibly after cascading), the total injection 
is 



Numerically, the dominant contribution to the injected photon energy comes from the most 
recent era, t > t eq , leading to a non-zero value for the Compton y parameter [100, 102J which 
quantifies deviations away from the black body spectrum. It is expected that the constraints 
on photon injection will be improved in the future by the ARCADE experiment [103] . 

Decaying cosmic string loops can also generate cosmic rays. The corresponding energy 
spectrum depends on the energies of the particles emitted in the loop decays. Recall that 
the fields making up the flat-direction strings consist of a light chiral supermultiplet and 
a heavy massive vector supermultiplet. In each cusp annihilation, both the heavy and the 
light states can be produced. The decays of the heavy states, with masses on the order 
of gv, can generate ultra- high-energy cosmic rays (UHECR) [104J. Decays of the light 
states, with masses on the order of m C ?;, contribute to the extragalactic diffuse gamma- 
ray background (EDGRB) |105j . To determine the relevant bounds and prospects, we will 
assume that the energy released in each cusp annihilation goes initially into a fraction Fi of 
the light states (with soft energy) and a fraction F h to the heavy states. We expect Fi ~ 1, 
with Fh possibly smaller. 

The contribution of decaying string loops to the EDGRB was studied in Ref. |105j . Data 
from the EGRET experiment [106] constrains the rate of energy emission into the light scalar 
states with masses on the order of 1000 GeV (that decay into lower-energy gamma rays) at 
the present time to dpi oop /dt Q < 4.5 x 10~ 23 eVcm -3 s _1 = 2.3 x 10~ 97 GeV 5 . Equating this 

14 The bound is fairly independent of the mass of the decaying particle. 
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bound with Eq. (15.161) evaluated at the present time, we obtain the bound [105 




2 x 10" 11 



Gn 



) 



< 1. 



(5.27) 



This does not represent a significant constraint beyond those found above. The heavy 
component states making up the string can also contribute to the EDGRB through the 
photons they produce in cascade decays. The limit in this case is about the same as from 
the decays of the light states given in Eq. fl5.27[ ). but with Fi replaced by F h . These constraints 
from the gamma-ray background on decaying cosmic string loops will be strengthened by the 
upcoming GLAST experiment |107j . However, the range of the model that can be probed 
may ultimately be limited by astrophysical background contributions to the gamma-ray flux. 

Ultra-high-energy cosmic rays can be produced by cusp annihilation if some of the heavier 
states making up the string are created. When the heavy states decay, their products are 
highly energetic, making them a source of high-energy neutrinos and UHECRs. Estimates of 
the UHECR flux for strings that decay into particles were made in Ref. [108] and are directly 
applicable to flat-direction cosmic strings. These authors find that for energies greater than 
about 6 x 10 9 GeV, the only relevant cosmic ray flux consists of neutrinos. The fluxes of 
highly energetic protons and photons are very suppressed because they are attenuated by 
their interactions with the cosmic background. Extrapolating the predictions of Ref. [108J, 
the neutrino signal from decaying strings can be probed directly at Ice Cube [109j down to 
G\i < 10~ 12 /F h in the energy range 10 5 GeV-10 8 GeV. The Auger project is sensitive to 
UHECR showers induced by energetic neutrinos in the energy range 10 9 GeV-10 11 GeV jllOj. 
The Auger measurements imply the constraint 



For Fh = 1, this corresponds to v < 10 13 GeV. 

Our analysis indicates that the visible matter signatures from decays of flat-direction 
string loops are consistent with observations provided Gfi is small enough. However, there is 
another visible matter signature that is challenging to reproduce in models of flat-direction 
strings, namely the baryon asymmetry of the universe. Flat-direction strings are formed 
following a period of thermal inflation. The typically low reheating temperature after 
thermal inflation, Eq. (14.61) . combined with the large amount of dilution from the inflationary 
expansion and reheating imply that baryogenesis mechanisms that operate at or above the 
electroweak scale will no longer work. Instead, the baryon asymmetry must be produced 
at very late times. This can arise from the strings themselves [1111 11121 I113j . from the 
non-thermal production of particles during reheating that have baryon-number violating 
decays [571 HIH HEl EDS, [111], or by the Affleck-Dine mechanism [TT8] . 

5.5 String Loops and Zero Modes 

In our discussion of radiation from cosmic string loops, we implicitly assumed that there do 
not exist any zero mode excitations along the strings. Zero modes are fermionic or bosonic 
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field fluctuations with vanishing energy that are localized on the string. The existence of 
zero modes on cosmic strings can alter the picture of loop radiation in important ways [119, 
IT20lfT2TlfT22] . 

These undamped, particle-like excitations can be excited when a string loop is formed. 
As the loop radiates and shrinks, the number density of the zero modes builds up. Eventually 
the angular momentum of the zero modes balances the tendency of the loop to shrink, and 
a quasi-stable loop remnant, or vorton, is left over. If such vortons are sufficiently long-lived 
and numerous, they behave like quasi-stable matter and can further modify the predictions 
of BBN or create too much dark matter. The presence of vortons typically leads to extremely 
strong constraints on the underlying field theory [122J. 

For the flat-direction strings we are studying, fermionic zero modes |123[ 1119] are of 
particular relevance. It was shown in Ref. [211 122] that such modes are a generic feature of 
super symmetric cosmic string solutions. In the present case, we also have supersymmetry 
breaking operators present in the Lagrangian. We find that adding a supersymmetry break- 
ing gaugino mass destroys all the fermionic zero modes. A recent study also suggests that 
more generally, fermionic zero modes do not form on closed string loops at all [124] . The 
existence of bosonic zero modes depends on the other fields in the theory and their couplings, 
and are less generic |119j . We do not consider them here. 

Zero modes, either bosonic or fermionic, are also unlikely to stabilize flat-direction cosmic 
strings simply because these strings are relatively wide. For the phase transition leading 
to flat-direction cosmic strings, we expect that the radius at which the zero modes would 
stabilize a string loop, if they were to exist, is usually much smaller than the width of the 
string [19]. As discussed above for cusp annihilation, when the separation between a pair of 
antiparallel string segments approaches the string width these segments will annihilate into 
their constituent fields, and the loop will decay before stabilizing as a vorton. 



5.6 Lensing by Cosmic Strings 

While an indirect gamma ray or gravitational wave signal from cosmic strings would be 
exciting, ideally one would like a direct observation to confirm their existence. This can 
be achieved by observing gravitational lensing by a string. The primary gravitational effect 
of the large mass density contained within a cosmic string is to modify the surrounding 
spacetime such that it is flat, but with a deficit angle of A6 = 87r G[i [125] . When light from 
a galaxy passes by a (non-relativistic) cosmic string, the deficit angle produces a distinctive 
double image with an angular separation of [31 1125] 

Aa = 8tt Gfi |p sin 0, (5.29) 

J-^os 

where Di s is the distance from the lensing string to the source galaxy, D os is the distance 
from the observer to the source, and sin <p is the angle between the string axis and the line- 
of-sight. From a single lensing event it is possible to determine Aa directly, as well as D os 
by measuring the redshift of the source. Given that a single string lensing event is found, it 
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is likely that the same string will also lens the images of other galaxies that are nearby on 
the sky [126j . By observing several galaxies lensed by the same string, the tension of that 
string can be determined [127] . 

The gravitational lensing signatures from flat-direction strings are even richer than those 
of ordinary strings because of the stability of higher winding modes. If many lensed images 
from different strings are observed, it may be possible to measure tensions of several strings 
and obtain clues about the mass spectrum of the higher winding modes. In this respect, flat- 
direction cosmic strings are similar to (p, q) cosmic superstrings. Both types of cosmic strings 
also have junctions connecting different winding modes. These can produce triple images, in 
addition to the double images produced by a lone string [128j . Since the spectrum of tensions 
of flat-direction strings is very different from that of (p, q) strings, the observation of many 
gravitational lensing events might allow one to distinguish between them. Unfortunately, 
the probability of observing a lensed image from a flat-direction string is very small due the 
indirect bounds on the tension, Gp < 10 -11 . This is smaller than the expected sensitivity of 
Gp ~ 10~ 8 from upcoming optical surveys [129j . and Gp ~ from the SKA [130] . radio 
survey [131j . The results of Ref. [127] also suggest that it will be difficult to determine the 
string tension accurately enough to distinguish between flat-direction string winding modes 
with similar values of the winding number N. 

6 Conclusion 

We conclude by summarizing some of our findings: 

• Abelian gauge symmetry breaking along a flat direction can give rise to strongly Type-I 
cosmic strings with tension p ~ 0.i7rt> 2 , gauge profile width of v~ l and scalar profile 
width w ~ m -1 , where mCt; characterizes the flatness of the potential. These flat- 
direction strings are likely to be formed after thermal inflation through flux-trapping. 

• The tension of the strings increases very slowly with their winding number N. Thus, 
higher-winding mode strings N = 2,3, . . . are energetically stable. This enables strings 
to be attracted to one another and zipper, creating stable formations with winding 
number N± + N 2 or \Ni — N 2 \, where N% and N 2 are the original string winding numbers. 

• Zippering affects the evolution of the resulting string network. Applying a simple 
network evolution model to flat-direction strings suggests that a large number of string 
modes develop roughly equal densities in the early universe. The total energy density 
is about the same as for a single string, but it is distributed among many species. 

• Flat-direction strings radiate gravitationally. However, in contrast to ordinary cosmic 
strings, they also may be able to radiate copiously into matter. The strings are expected 
to fully radiate away, as there is no vorton obstruction for the supersymmetric flat- 
direction strings under consideration. 
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• In contrast to GUT strings, flat-direction strings are generically compatible with cur- 
rent direct observational constraints, Gji < 3 x 10~ 7 [63j [MJ ESI [661 [67] . If the typical 
initial loop length is close to the horizon scale, LISA and upcoming millipulsar timing 
probes may be able to detect the gravitational wave signal from these strings. However, 
the gravity wave signal at higher frequencies is suppressed for flat-direction strings, 
making their detection at LIGO extremely challenging. 

• Particle emission from cusp annihilation is likely to be the dominant loop decay 
mechanism if the loop length is always much smaller than the horizon. This intriguing 
prospect for flat-direction cosmic strings entertains the possibility that ultra-high- 
energy cosmic rays or nonthermal dark matter originate from their particle emission. If 
all loops decay entirely into particles, the constraints from BBN, the CMB blackbody, 
and UHECRs imply the bound v < 10 13 GeV, corresponding to G/i < 1CT 13 . 

• The multi-tension network of flat-direction strings formed in the early universe is in 
contrast to the standard single-tension string networks, but similar to (p, q) cosmic 
superstring networks, and thus may mimic the latter by giving rise to multiple lensing 
events. However, the spectrum of tensions of flat-direction strings is constrained by 
indirect bounds, and may be too low to be observed in the near future. 

We find that flat-direction cosmic strings behave in ways that are qualitatively different 
from both ordinary (abelian Higgs) cosmic strings as well as (p, q) cosmic superstrings. These 
differences in behavior may be distinguishable through probes of the early universe. 
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